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Abstract 

We introduce a global scheme on the n-torus of a controlled incom- 
pressible Navier-Stokes equation in terms of a coupled controlled infi- 
nite ODE-system of Fourier-modes with smooth data. We construct 
a scheme of global approximations related to linear partial intcgro- 
diffcrcntial equations in dual space which are uniformly bounded in 
dual Sobolev spaces with polynomially decaying modes. The scheme 
is based on some infinite matrix algebra related to weakly singular in- 
tegrals, and a Trotter-product formula for dissipative operators which 
leads to rigorous existence results and a uniform bound for the solutions 
of the successive approximating linearized equations in dual space. For 
data with polynomial decay of the modes (equivalent to smooth data 
on the n-torus) global existence follows from a compactness property 
in strong dual Sobolev spaces (spaces of polynomially decaying Fourier 
modes). The difference to schemes considered in [5J [7J EH E] is that 
we have a priori estimates for solutions of the iterated global linear par- 
tial intcgro-diffcrential equations. Furthermore, the control function 
does only control the non-dissipative zero modes and is, hence, much 
simpler. The analysis of the scheme leads to an algorithmic scheme of 
the Navier-Stokes equation on the torus for regular data (certain poly- 
nomially decaying modes) based on the mentioned Trotter-type prod- 
uct formula for specific dissipative operators. This algorithm which 
is closely related to the Trotter product formula is natural, since the 
latter is a limit formula with respect to the order of the modes and 
expresses an approximation of approximating functions via finite sys- 
tems. It is our main intention here to define an algorithm wich uses the 
damping via dissipative modes and which converges in strong norms 
for arbitrary viscosity constants v > (or Reynolds numbers). Finally, 
we discuss some notions of the concept of 'turbulence' which is not an 
effect of a blow-up of vorticity, but may be a concept to be described 
by dynamical properties, as has been suggested by other authors for a 
long time (cf. 013], and @]). 

2000 Mathematics Subject Classification. 35K40, 35Q30. 
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1 Introduction 



The formal transformation of partial differential equations to infinite sys- 
tems of ordinary differential equations via representations in a Fourier basis 
may be useful if the problem is posed on a torus. Formal representations of 
solutions are simplified at least if the equations are linear. Well, even non- 
linear infinite ODE-representations may be useful in order to define solution 
schemes. However, in order to make sense of an exponential function of 
an infinite matrix applied to an infinite vector the sequence spaces involved 
have to be measured in rather strong norms. Schur [TB] may have been the 
first who provided necessary and sufficient conditions such that infinite lin- 
ear transformations make sense. His criteria were in the sense of ^-sequence 
spaces and are far to weak for our purposes here. However, dealing with 
the incompressible Navier-Stokes equation we shall assume that the initial 
vector- valued data function h = (hi, • • • , h n ) T satisfies 

hi G C°° (T n ) , (1) 

where T n denotes the n-torus (maybe the only flat compact manfold of 
practical interest when studying the Navier-Stokes equation). In the fol- 
lowing let Z be the set of integers, R the field of real numbers, and let 
a, (3 G Z n be multiindices. Differences of multiindices a — (3 G Z n are built 
componentwise, and as usual for a multiindex 7 with nonnegative entries 
we denote | -y ] = ^27=1 17*1- Consider a smooth function h = (hi, • • • , h n ) T 
with hi G C°° (Tf), where Tf = W 1 jVU 1 is the torus of dimension n and size 
I > 0. Note that the smoothness of the functions hi means that there is a 
polynomial decay of the Fourier modes, i.e., for a torus of size I = I for the 
modes of multivariate differentials of the function hi we have 

-D 7 /ii Q := J Jn D%hi(x) exp (— 2-niax) dx 

= (-l)\^\ J Jn /i i (x)n™ =1 (-2vria J ) 71 exp (2iriax) dx (2) 

= U^ =1 (-2Triai)^h ia , 

where as usual for a multiindex 7 with nonnegative entries 7i > the symbol 
Dx denotes the multivariate derivative with respect to x and of order 7$ 
with respect to the ith component of x. Here hi a is the ath Fourier mode 
of the function hi, and D^hi a is the ath Fourier mode of the function Djhi. 
Since Djhi is smooth this function has a Fourier decomposition on T n . The 
sequence (-D 7 ^i Q ) Q , gZ n of the associated a-modes exists in the space of square 
integrable sequences I 2 (Z n ). Hence, we have polynomial decay of the modes 
hi a as M t 00 • This polynomial decay is important in the following in order 
to make sense of certain matrix operations with certain multiindexed infinite 
matrices and multiplications of these matrices with infinite vectors, because 
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we are interested in regular solutions. In the following we have in mind that 
a multiindexed vector 

u F := (u a )l €Zn (3) 

(the superscript T meaning 'transposed') with complex constants u a (al- 
though we consider only real solutions in this paper) which decay fast enough 
as | a | f 00 corresponds to a function 

u G C°° (Tf) , (4) 

where 

u(x) := ^ u a exp ( — - — j. (5) 

If the representation can be rewritten with real functions cos ( ^lax ^ an( j 
sin ( 2 ™ x ) and with real coefficients of these functions, then the function in 
@ has values in R, and this is the case we have in mind in this paper. We 
say that the infinite vector is in the dual Sobolev space of order s £ R, 
i.e., 

u F G h s (II 1 ) (6) 

in symbols, if the corresponding function u defined in @ is in the Sobolev 
space H s (Tf)of order s£R. We may also define the dual space h s directly 
defining 

u F G h s (Z n ) ^ l u «| 2 < °°> ( 7 ) 

where 

(a) := (l + | a | 2 ) 1/2 . (8) 

The two definitions are clearly equivalent. 

Note that being an element of Q for nonnegative s means that we 
have some decay of the a-modes. For example for constant a-modes the 
corresponding function in classical -£P-space 

E Cexp(^) (9) 
is a formal expression of C times the 5 distribution and we have 

Ti 

5 G H s if s < --. (10) 

Note that we may consider n-tori which are build by identification of oppo- 
site sides of the n-cube [0, l] n not [—1, l] n . This means that we have Fourier 
representations in the space of even functions, i.e., we may assume that 
u_ a = u a . We make this remark because it is a useful fact for computa- 
tion, but it is not needed for the proof of our main theorem. Now we may 
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rephrase the ideas in [7], [8], |10j . and |1 1 j in this context. We shall devi- 
ate from these schemes as we simplify the control function and define an 
iterated scheme of global equations corresponding to linear partial integro 
differential equations approximating the (controlled) incompressible Navier 
Stokes equation. Formally, the modes (vi a ) a ^z n , 1 < i < n, oi the velocity 
function Vi, 1 < i < n of the incompresssible Navier-Stokes equation satisfy 
the infinite ODE-system (derivation below) 



dv ia _ v^n [ _47ro| \ _ „ n „ 

(11) 

where the modes V{ a depend on time t and such that for all 1 < i < n and 
all a 6 Z n we have 

v ia {0) = h ia . (12) 

In the infinite system and for fixed a we call the equation in with left 
side ^f - the a-mode equation. Some simple but important observations are 
in order here. First note that the damping 




47ra] 



(13) 



is not equal to zero unless \a\ = Yl7=i l a *l = ^- Hence we have damping 
except for the zero-modes Vio (where the subscript denotes the n-tuple 
of zeros). Second, note that the zero modes Ujo contribute to the a-mode 
equation for a / only via the convection terms. This is because for the 
second term on the right side of (|lip only the term with 7 = a, i.e., the 
summand 

l-^-VjoVia (14) 

5=1 



contains a zero mode. Third, note that the pressure term in (|lip has no 
zero mode summands since 



2mail {am YJU^ = 

2ma i k»m V- : w 



(15) 



Furthermore, note that the 0-mode equation consists only of terms corre- 
sponding to convection terms, i.e., we have 

= ~ Y?j=\ E 7 gZ"\{o} ~T' v j(--/) v n- ( 16 ) 
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Furthermore we observe that in (|16p we have no zero modes on the right 
side but only non zero- modes of coupled equations. Furthermore all non- 

zero modes involve a damping term, because v I jr- > for v > 



IF 

and \a\ 7^ 0. Note that damping becomes stronger as the order of the 
modes \a\ 7^ increases - an important difference to the incompressible 
Euler equation, where we have no damping indeed. This is also a first 
hint that it may useful to define a controlled incompressible Navier-Stokes 
equaion on the torus where a simple control functions controls just the zero 
modes. Next we define schemes for computing the modes. For purposes of 
global existence it can make sense to start with the the multivariate Burgers 
equation, because we know that we have a unique global regular solution for 
this equation on the n-torus. This can be derived by the a priori estimates 

^\Ht,.)\\ H s<\\u(t,.)\\ H3+ iJ2 J2 \\D a u i D^u j \\ L 2-2\\Vu\\ 2 HS , (17) 

i,j \a\+\P\<s 

which hold on the n-torus, or by the arguments which we discussed in |9] 
and [TU]. In our context this means that for positive viscosity v > and 
smooth data (at least formally) we have a global regular solution for the 
system 

Fit = YTj=i v ( IF' J Uia ~ S™=i E 7 ez« ^FF Ku j{a--y) u i'y ( 18 ) 

where V{ a depend on time t and such that for all 1 < i < n and all a£Z™ 
we have 

u ia (0) = h ia . (19) 

This suggests the following first approach concerning a (formal) scheme for 
the purpose of global existence. We compute first v® a = ui a for 1 < i < n 
and a£Z n and then iteratively 

-l+J2 Sv i*> ( 2 °) 
P =i 



where 5vf a := v v iOL — i?- 1 for all 1 < i < n and a£Z" and for p > 1 we have 

(21) 

+2mal 



En 1 j \ P v-mi v-v 



E"fc=lE 76 zn 47r 7j(«fc-7fe)«f 1 



where i>i a depend on time t and such that for all 1 < % < n and all a G I/ 1 
we have 

«ia(0) = h ia . (22) 
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Alternatively, we may start with the initial data hi and their modes hi a as 
first order coefficients of the first approximating equation. In this case we 
have for an iteration number p > the approximation at the pth. stage via 
the linear equation 



dt — 2^j=l u [ -1^ V ia Ljj=\ Z^-yeZ n l U j (0-7)^17 

(23) 



+2mai l{ a ^o} i^f ' 

where again for all 1 < % < n and all a £ Z" we have 

^«(0) = /w (24) 

For p = we then have \ = ^ := /i J '( Q ,_ 7 ). The latter equation 

in (|23p still corresponds to a partial integro-differential equation. However, 
this scheme has the advantage that we can built an algorithm on it via a 
Trotter product formula for infinite matrices (because we know the data hi a 
and for p > 1 the data vf^ 1 from the previous iteration step. 

Remark 1.1. The equation in (|2ip corresponds to a linear integro-differential 
equation in classical space and differs in this respect from the approximations 
we considered in [10] and [11] and also in [8] and [7] (although we mentioned 
this type of global scheme in [7]). The analysis in all these papers was 
based on local equations because a priori estimates are easier at hand - as 
is the trick with the adjoint. In dual spaces of Fourier basis representation 
considered in this paper spatially global equations can be analysed easier. 
We could also use the corresponding local equations. The corresponding 
system is 



. o • -, Ei,fc=iE 7e zn47r 7i (a fe - 7fe )^ 7 1 ^ 1 

+27^1(^0} TJU^k 



(25) 



but there is no real advantage analyzing (|25|) instead of ([211 , 
Note that the 'global' term in OTj) is 

. . , Ei,fe=i E 7G Z" 47r 7j(afc - 7fc)^j 7 " 1?; fc(a-7) , 0(! x 
27rzaa { ^o} w , (26) 

and it looks more like its 'local' companion in this discrete Fourier-based rep- 
resentation than is the case for representations in classical spaces. Formally, 
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for all 1 < i < n and q£Z" and for p > 1 for 5vf a = vf a — vf a 1 , a £ Z n we 



have 



En / 47rQ j \ r p v^n v-^ 2iri"/j p—1 p 

i=l V I W J 0V ia ~ 2^j=l Z^ 7 GZ« I V j( a -~ ( ) ov i'Y 



(27) 

+2mai 1 W0 } E? w 

££fc=i E 7S Z" 47r 7j(«fc-7fc)^ 2 <^ 1 _ 7) 



-2mail {am „„ 2 



where Vj a depend on time i and such that for all 1 < i < n and all a 6 Z n 
we have 

fo; Q (0) = 0. (28) 

It is well-known that the viscosity parameter may be chosen arbitrarily. This 
fact facilitates the analysis a bit, but for the proof of global existence it is 
not essential. It is important that the viscosity is strictly positive [y > 0). 
We need this for the Trotter product formula below, and without it the 
arguments breaks down. It is also useful to have a stronger damping for 
simulations- so let us be more explicit about this. Consider a parameter 
transformation of the form 

Vi (t,x) = v*(r,y) = r x v*(r%r»x), p(t,x) = p*(r,y) = r s v* (r"' t,r»x) 

(29) 

for some positive real number r > and some positive real parameters 
X,fj,,v'. We have (using Einstein notation for spatial variables) 

&Vi(t,x) = r x +»'£v*(T,y), v id (t,x) = r A+ ^(r,y), 

(30) 

Vi,j,k{^ x ) = rA+2M &^T<( r >2/)' P,i& x ) = r 5+ Vp%T,y). 

Hence 



" X+ " Q^v*(T,y) = vr x+2 »Av*(T,y) -Y,r 2X+ »v*(T,y)^(T,y) =r S ^ 

(31) 



j= i - 3 



which (for example) for the parameter constellations 

fi, i/, A, 5 with A + fi — v = and 5 + {/,— A — z/ = (32) 
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becomes 

ur 2 »-»'Av*(T, y) - £™ =1 r x +^'v*(r, y)^(r, y) - r^'^'p^ (33) 



z,r 2 ^'A<(r,y) - £? =1 «*( r ,y)§|(T,y) 



P.. 



Note that we may choose (for example) 2[i — v' G K + (]R + being the set 
of strictly positive real numbers) freely and still satisfy the conditions (|32|) . 
Note that we can take v' = v. Hence, we are free to assume any viscos- 
ity v > 0. Concerning effciency of the corresponding algorithm such an 
observation can be useful, although it is of limited use when it comes to 
more realistic models of fluids with boundaries and free boundaries. For the 
numerical and computational analysis it is useful to refine this a bit and 
consider certain bi-parameter transformation with respect to the viscosity 
constant v > and with respect to the diameter I of the n-torus. From an 
analytic perspective we may say that it is sufficient to prove global existence 
for specific parameters v > and / > 0. Let us consider coordinate trans- 
formations with parameter / > first. If v \p l is solution on the domain 
[0,oo) x T™, then the function pair (t, x) — > v l (t,x), and (t,x) — > p l (t, x) on 
[0, 00) x T" along with 

v\t,x):=\^(^ |V (34) 

and 



p\t,x):=^p l (^-^j (35) 
is a solution pair on the n-torus of size I > with initial data 

h l (x):= yhHfjonT?. (36) 

Hence if we can construct solutions to the Navier-Stokes solution for hi G 
C°° (T™) without further restrictions on the data h;, then we can construct 
global solutions for the problem in with h = hi G C°° G (T n ) 1 with 
factor / as in (|36p . We only need to produce a fixed number I > such that 
arbitrary data are allowed. Second if v\ , 1 < i < n is a solution to the 
incompressible Navier Stokes equation with parameter v > then via the 
time transformation 

Vi(t,x) := (rv)~ 1 v^ ((ru)~ 1 t,x), p(t,x) := (ri/) _2 p I/ ((r^) _1 t, x) (37) 

(for all t > 0) we observe that v i is solution of the incompressible Navier- 
Stokes equation with viscosity parameter which may be chosen. Concatena- 
tion of the previous transformations shows that we can indeed choose specific 
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values for v > and I > 0, and this may be useful for designing algorithms. 
For example, it is useful to observe that in (|27p the modulus of the diagonal 
coefficients 



St^J (38) 



becomes large for \a\ 7^ if v is large or I is small in comparison to the other 
terms of the iteration in (|27|) . but it is even more intersting that we may 
choose v large and I > large on a scale such that the convection terms are 
small in comparison to the diagonal damping terms in ([38|) . As we observed 
in other articles for the time-local solution of the Navier-Stokes equation we 
may set up a scheme involving simple scalar linear parabolic equations of 
the form ([39]) (cf. [9l QUI [JT] ) . We maintained that in classical space a 
global controlled scheme for an equivalent equation can be obtained if in a 
time-local scheme the approximate functions in a functional series (evaluated 
at arbitrary time t) inherit the property of being in C k n H k for k > 2 
together with polynomial decay of order k. This is remarkable since you 
will not expect this from standard local a priori estimates of the Gaussian. 
For the equivalent controlled scheme in [11] we only needed to prove this 
for the higher order correction terms of the time-local scheme. There it is 
true for some order of decay because the representations involve convolution 
integrals with products of approximative value functions where each factor 
is inductively of polynomial decay. The growth of the Leray projection term 
for the scheme in [11] is linearly bounded on a certain time scale. In [7J[8] we 
proposed more complicated equivalent controlled schemes with a bounded 
regular control function which controls the growth of the controlled solution 
function such that it stays bounded for all time. Actually, the simple schemes 
may be reconsidered in terms of the Trotter product formula representations 
of solutions to scalar parabolic equations of the form 

— -vAu + Wu = g, (39) 

with initial data u(0,x) = f 6 n sg KnfP and some dynamically gener- 
ated source terms g, and where W = Y^a=i w i( x )~Bx~ ^ s a vec t° r held whith 
bounded and uniformly Lipschitz continuous coefficients. The Wi then are 
replaced by approximative value function components at each iteration step 
of course. Local application of the Trotter product formula leads to repre- 
sentations of the from 



exp (t (uA + W))f= lim j^exp I - W J exp (^ A J J /■ ( 40 ) 

Then in a second step having obtained time-local representations of solutions 
to the incompressible Navier-Stokes equation one can use this reduction from 
a system to a scalar level and apply the Trotter product formula on a global 
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time scale. Prom this point of view the local contraction results considered 
in [TJ] or in [9] (with the use of adjoint equations) are essential. In dual 
spaces the use of a Trotter product formula seems to be not only useful but 
mandatory in two respects. Consider the infinite matrix 



D ■■= ( d *ft)«,B&» ■= <w y; v\ — (4i) 




with the infinite Kronecker delta function 5 a p, i.e, 

1 if a = p 



(42) 



if a ^ 0. 



If we measure regularity with respect to the the degree of decay of entries 
as the order of the modes increases, then this matrix lives in a space of rather 
weak regularity. Worse then this, iterations of the matrix, i.e., matrices of 
the form 



(43) 



live in matrix spaces of lower and lower regularity as m increases (if we 
measure regularity in relation to decay with respect to the order of modes). 
However, since the matrix D has a minus sign we may make sense of the 
matrix 



exp(D) = (6 a p exp (d aa )) a &&n = 8 a pexp 




(44) 

and this matrix makes perfect sense in terms of the type of regularity men- 
tioned, i.e. the type of regularity expressed by dual Sobolev spaces which 
measure the order of polynomial decay of the modes. Since we consider 
infinite ODEs involving a Leray projection term we shall have correlations 
between the components 1 < i < n even at each approximation step. Hence 
our matrices shall involve big diagonal matrices of the form 

(45) 

where tfy denotes the usual Kronecker 5 for 1 < i,j < n and exp(D) is as 
in (|44|) above. However, the idea that such a dissipative matrix lives in a 
regular (='polynomially decaying modes as the order of modes increases') 
infinite matrix space is the same. 
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Now consider the first step in our scheme for the vectors (vf a ) aeZ „- For 
p = the modes {vi a ) a&Zn are equal to the modes of the controlled scheme 
we describe below and the corresponding equation can be written in the 
form 



dt 



(46) 



O • 1 E",fc=l E 76Z " 47r 7j(Qfc-7fc)^° 7 fefc(c t - 7 ) 



If we consider the n-tuple v F = (v F , • • • , v^) T as an infinite vector with 
vf := (via) aeZ n, then with the usual identifications the equation ((15]) is 
equivalent to an infinite linear ODE 

V-^- (47) 

where the matrix is defined implicitly given by (|46|) and will be given 
explicitly in our more detailed description below. Together with the initial 
data 

v°' F (0) = h F (0) (48) 

this is an equivalent formulation of the equation for the first step of our 
scheme. We shall define a dissipative diagonal matrix Dq and a matrix Bq 
related to the convection and Leray projection terms such that Aq = Dq+Bq, 
and prove a Trotter product formula which allows us to make sense of the 
formal solution 

v°' F (t) = exp(A t)h F 

(49) 

:= lim itoo lim fctoo (exp (P M i( D o)i) exp (P M iB Q ±)) h F . 

Here P M i denotes a projection to the finite modes of order less or equal to 
I > 0. Here and in the following we may assume that we have some order of 
the multiindeices ant that this order is preserved by the projection operators. 
Note that at this first stage the modes h{ a are not time-dependent. Hence 
Aq is defined as a matrix which is independent of time, and we have no need 
of a Dyson formalism at this stage. For the higher stages of approximation 
a control function comes in and we have to deal with time dependence of 
the related infinite matrices A p which define the infinite ODEs at iteration 
step p > for the modes vf a and v r ^ in the presence of a control function 
r. The formula ()40|) makes clear why strong contraction estimates of the 
time-local expressions 

exp Uw\ exp (J-uA\ f (50) 
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are important. As we said, we pointed out this in [10] and [TT] from a differ- 
ent point of view. In dual spaces we may approximate such expressions via 
matrix equations of finite modes (projections of expressions of the form (|50p 
to approximating equations of finite modes) . For the finite mode approxima- 
tion of the first approximation sequence v® a , a € Z n first order coefficients 
are time independent and we may use the Baker-Campbell-Hausdorff for- 
mula for finite matrices A and B of the form 

exp(A) exp(B) = exp(C), (51) 

where 

C = A + B + ±[A,B] + ^[A, [A,B]] + ±[[A,B],B] 

(52) 

+ higher order terms. 

Here [., .] denote Lie brackets and the expression 'higher order terms' refers 
to all terms of multiple commutators of A and B. We shall see that for dissi- 
pative operators as in the case of the incompressible Navier-Stokes equation 
and its linear approximations we can prove extensions of the formula in (|5ip 
to infinite matrices, where we restrict our investigation to the special cases 
which fit for the analysis of some infinite linear ODEs approximating the 
incompressible Navier-Stokes equation written in dual space. Note that for 
higher order approximations vf a , a£Z", 1 < i < n with p > 1 we have time 
dependence of the coefficients and this means that we have to apply a time 
order operator as in Dyson's formalism in order to solve the related linear 
approximating equations formally. Note that (|52|) is closely connected to 
the Hormander condition, and this was one of the indicators which leads to 
the expectation in [8] that global smooth existence is true if the Hormander 
condition is satisfied. This will be considered elsewhere. 

In this paper we shall see that we can simplify the schemes formulated 
in classical spaces in our formulation on dual spaces in some respects. The 
first simplification is that we may define an iteration scheme on a global time 
scale. The second simplification is that the estimates in dual spaces become 
estimates of discrete infinite sums which can be done on a very elementary 
level. On the other hand we shall need a control function which ensures 
that the scheme for the controlled equation is a scheme for non-zero modes 
(ensuring that the damping factors associated with strictly positive viscos- 
ity or dissipative effects are active for all modes of the controlled scheme). 
We shall choose specific v and I and define a controlled scheme for a con- 
trolled incompressible Navier-Stokes equation such that the diagonal matrix 
elements corresponding to the Laplacian terms become dominant and such 
that there exists a global iteration in an appropriate function space. How- 
ever the control function is a scalar univariate function and is much simpler 
than the control functions considered in |7] and [8]. Indeed, the control func- 
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tion in the present paper will only control the zero modes Vio of the value 
function modes. 

Next we define the controlled scheme. We may start with one of the two 
possibilities mentioned above. We may start with the solution scheme for 
the multivariate Burgers equation, which is given in dual representation by 
the infinite ODE 



du ia _ v^n. I _4tto| \ _ „ n „ /ro\ 

dt ~ 2^j=l u \ —p—jUia 2sj=l Z^7GZ™ I u j(a--y) u i-y> K 06 ) 



where the modes U{ a depend on time t and such that for all 1 < i < n and 
all a G Z n we have 

u ia (0) = h ia . (54) 

This is a slight variation of the linearized scheme mentioned above. We may 
define vf a = Ui a for 1 < i < n and a G Z n and then iteratively 

k 

P =i 

where for p > 1 we define <5t>f Q := ff a — uj^ ^ for all 1 < i < n and q£Z" 
such that 



dt - 2_,j=l V I j5- I - 2^=1 2^ 7 eZ"\{a} I V j( a - 1 ) 0V i-i 

Ei=i ^FKo" 1 ^" 1 (56) 



+ 1 {q ^ } E n =i 4TO 2 



. „ . -. E"fc=iE 7eZ «47r 7j (a fc -7 fe )5^ 7 1 ^ (Q 1 _ 7) 
+^ 7rta i i {«5«)} EIU4^ 

(where for p = and p = —1 we define 5vf a = 0). Note that we extracted the 
zero modes from the sums in (|56p . The reason is that we want to extend this 
scheme to a controlled scheme which is equivalent but has no zero modes. 
The present local scheme suggests this because the Leray projection terms 
do not contain zero modes. Note that for all 1 < i < n and all a 6 Z" we 
have 

5*4(0) =0. (57) 

Global smooth existence means that the vf a are defined for all time = 
{t G R|t > 0}, and such that polynomial decay of the modes is preserved 
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throughout time. Next we introduce the idea of a simplified control function 
was introduced in [7] and [8], and [11] in a simplified but still complicated 
form. Here we introduce a control function for the zero modes. We define 

< a P (t) = vl(t) + r p (t), (58) 

where 

r p : [0, oo) E (59) 

is defined by 

r P (t) = -«£(<)- (60) 

We have to show then that ro : [0, oo) — )■ R is well-defined (especially 
bounded), i.e., we have to show that there is a limit 

oo 

r (i) := hmrg(t) = rg(t) + £>r?(i) (61) 

p=l 

along with ^^(t) = r^i) — rQ _1 (t). This leads to the following scheme. We 
start with the solution for the multivariate Burgers equation but annihilate 
the zero modes, i.e. we define 

v ia = v ia = u i<* ( 62 ) 

for a/0 and 

«a° = «a + = o, (63) 

where for all t > 

rg(t) = -««)(*). (64) 
For p > 1 and for 1 < i < n and a £ 7* n \ {0} we define 

r '° + (65) 
p=i 



v ta 



where £ujf := + r p a for all 1 < i < n and a£Z"\ {0}, and 

_ V n u ( ^ a J I Ai, r ' p V n V 27ri 7j r,p-l <- r,p 

i i (66) 

+27riail {Q#0 } £ 7r 74^ 

En v-^ /i / \ r r.p— 1 r.p — 2 

+^*Oi A {a#0} V" I W ' 

^11 — 1 £ 
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Furthermore, for all 1 < i < n and a = we shall ensure that 

6v%(0) = 0. (67) 
Note that in the equation (j66|) the terms 



" E"=l ^O 1 ^ - E •=! 1 (68) 

on the right side of (|56p are cancelled. This is because we define the control 
function r$ such that the zeromodes become zero. This is done as follows. 
First we note that for p = and p = —1 we may define Sv^ = for all 
a € Z n . For a = we define first the increment <5w*g r ' p for p > 1 via the 
equation 

dSv M T ' P _y" , . / f™l \ A *' r ' p V" V 27r» 7j *,r,|)-lr *,r,p 

_ v^n 2xvjj r *,r,p-l *,r,p-l 

Z^j = l 2^jeZ n \{a} I 0U j(-~ f ) u i-y 

(69) 

Then we define 

K(t) = -5v*™, (70) 
closing the recursion. Note that this ensures 

Hi = o (7i) 

at all stages p > 1 for the zero modes (provided that we can ensure that rg is 
globally well-defined) . The introduction of the control function for the zero 
modes makes the system 'autonomous' for the modes a 6 Z n \ {0}. This 
means that all modes have damping factors (involving v > 0), i.e., we have 
the stabilizing diagonal factor terms 



(72) 




which contribute for 

ol) / 0, or oi / V 1 < i < n. (73) 

3=1 

We observed that the parameter v > may be large without loss of gener- 
ality. This alone may indicate that the matter of global smooth existence 
should be detached from the subject of turbulence as simulation indicate 
turbulent phenomena for high Reynold numbers. Solutions of each approx- 
imation step of the scheme can be represented in terms of fundamental 
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solutions of scalar equations, where each iteration step of higher order in- 
volves fundamental solutions of linear scalar parabolic equations. This is 
an idea which we considered earlier in [9], [TO], [8], and [7]. In this paper 
we consider 'fundamental solutions' in a dual space of Fourier modes. They 
may be called 'generalized ^-functions' but in general they live only in dis- 
tributional space of negative Sobolev norm. However, we shall see that for 
certain initial data we can make sense of the related analytic vectors, and 
show that there spatial convergence radius is global (holds for the full size 
of the n-torus). 

In this paper we observe that for all data h = (hi, ■ ■ ■ , h n ) along with 
hi € C°° there exist a number v > such that the global extended iteration 
scheme converges globally to a solution of an infinite ODE system equiv- 
alent to a controlled Navier-Stokes equation in a strong norm, i.e. for all 
t > the sequences Vi a , q£Z™ converge in the dual Sobolev space h s (Z n ) 
for s > 2. Moreover, the control function is a globally well-defined uni- 
variate bounded differentiable function such that controlled incompressible 
Navier-Stokes equation is equivalent to the usual incompressible Navier- 
Stokes equation system. More precisely, we have 

Theorem 1.2. There are numbers I > and v > such that for all 1 < 
p<n and data h G [C°° (Tf )] n for alls G R we have 

< a (t)=lim<nt)e^(P). (74) 

mjoo 



for the scheme described in the introduction. Moreover, t — > v\ a (t) satis- 
fies the infinite ODE system in [56\) . This implies that Vi a (t) = v\ a (t) — 
b~0a r i{t) satisfies the infinite ODE-system corresponding to the incompress- 
ible Navier-Stokes equation, where 5o a = i/a / and 5 a Q = 1 if a = 0, 
i.e., the function 

Vj := Vj a exp (2iriax), 1 < j < n, (75) 

aeZ n 



is a global classical solution of the Navier stokes equation system J76|) below. 



Note that Vj a , a G Z" are functions of time as is vj in (75). 

Here we say that v = (v±, ■ ■ ■ ,v n ) T is a global classical solution of the 
incompressible Navier-Stokes equation system 

' |f - + (v • V)v = - Vp t > 0, 
V • v = 0, t > 0, , (76) 

v(0, x) = h(x), 

on some domain f2 (which is the n-torus T n in this paper) if v solves the 
equivalent Navier-Stokes equation system in its Leray projection form, where 
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p is eliminated by the Poisson equation 

n 

- Ap = v j:k v k ,j. (77) 
i,fc=i 

Note that Vj^ denotes the first partial derivative of Vj with respect to the 
component etc. (Einstein notation). This means that we construct a 
solution to an equation of the form 

( |-i/Av + (v- v)v = z]k=i I VK ( X - vKkVkA^ y) d y 

(78) 

^ v(0, x) = h(x), 

along with ([77]) . For a fixed function v which solves ([78]) with (fTTj) the 
Cauchy problem for divergence 

^divv - i/A div v + J2]=i v j div v + YTjk=i v j,kVk,j = - Ap 

(79) 

divu(0, .) = 
simplifies to 

Jjdiv v — vA div v + Y^J=i v j div v = 

(80) 

div v(0, .) = 0, 

yielding divv = as the unique solution of (|80p . For this reason it suffices 
to solve the Navier-Stokes equation in its Leray projection form (as is well- 
known). In the next Section we look at the structure of the proof of the 
main theorems. Then in Section 3 we prove theorem 11.21 

Corollary 1.3. The statements of hold for all numbers I > and 

v > 0. Furthermore, the solution is unique in the function space 

F ns = {g : [0,oo) x W 1 -> R n \ 9i G C 1 ([0,^),H S (R n )) & 9i (0, .) = hi} , 

(81) 

where g = (g\, ■ ■ ■ , g n ) T and hi G H s for s > n > 3. 

Remark 1.4. For the assumptions of (jl.3p the methods of this paper lead 
directly to contraction results. Concerning uniqueness the regularity as- 
sumptions of the initial data may be weakened, of course. 

2 Structure and ideas of the proof of theorem 11.2 
and corollary 11.31 

In the proof of theorem 11.21 in section 3 below we first recall that the in- 
compressible Navier-Stokes equation is formally equivalent to a system of n 
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infinite ODE-systems 

dv ia _ V n „(_ A ™*\„. _ v n V 2 " 7 ^ - , 



(82) 



o • -i E"fc=i E 7 ez" 47r 7j("fe-7fe)^'7' ,; fc(c«-7) 

-27rzai l {a#0} — v? : W 



where a£Z", and where the initial data vf (0) = (vi a (0)) aeZn for 1 < i < n 
and a£Z° are given by n infinite vectors 

vf(0)=hf = (M T - (83) 



Furthermore, the entries /ij Q of the latter infinite vector denote the Fourier 
modes of the functions hi along with h £ [C°° (Tf)] n . Smoothness of h 
translates to polynomial decay of the modes hi a . Here we may say that the 
modes hi a have polynomial decay of order m > if 

\a\ m h ia J, as \a\ = V] \cti\ f oo (84) 

i=l 

for a fixed positive integer m, and the modes have polynomial decay if 
they have polynomial of any order m > 0. We may use the terms 'modes of 
solution have polynomial decay' and 'solution is smooth' interchangeably in 
our context. We may rewrite the equation (j82[) in the form 



*£ = A N S(v)v F , (85) 
where = (vf, • • • , v^) T . Furthermore A NS (v) is a nil 1 x raZ n -matrix 

^(v) = «(v)) 1 ^ n (86) 
where for 1 < i,j < n the entry Af- S (v) is a Z n x Z n -matrix. We define 

In n \ 

a ns (v) v F =lj2 A ™ (-) v f . • • • > £ ( v ) v f > ( 87 ) 

where for all 1 < i < n 

The entries A^^ (v) of (v) are determined by the equation in (|82p of 
course. On the diagonal, i.e., for i = j we have the entries 



<MS (v) = <% E?=i - (-^f) - s* E -=i 

Z-^i— 1 i 



n 27r^ 7 - 



(89) 
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where Sij = 1 iff i = j and zero else denotes the Kronecker <5-function, and 
off-diagonal we have for i ^ j the entries 

(1 - 5 lj )A iajp (v) = 27rza i l {Q ^ 0} ^ n 4 ^ q2 ■ ( 90 ) 

Here, we remark again that the i in the context 2ni denotes the complex 
number i = \J — 1, while the other ith are integer indices. The next step is 
to define a contolled system according to the ideas described in the intro- 
duction. Recall the main idea: the control function cancels the zero modes 
in an iterative scheme. If this is possible and a regular limit exists, then it 
cancels the zero modes in the limit as well. Consider the representation in 
(|82|) . The zero modes do not appear in the dissipative term related to the 
Laplacian, and we have observed that we may represent the Leray projection 
term without zero mode terms as well. Hence, for the zero modes viq in the 
equation in (j82|) we have for a = 

dt ~ 2sj=i Z^ez™ i u j(--y) u n- y UL > 
Assuming that the function t — > Vio(t) is bounded it is natural to define 

rg(t) := -v i0 (t). (92) 

Formally, this leads to the controlled equation for v\ a , a G Z n \ {0} = Z n '°. 



dt — 2-y=l v I ~T~ I V ia l_ij=\ 2^ 7 gZ"\{0,q} I v j( a --y) v ij 

+27rza i l{ c ^o} v" : W ' 



(93) 



where a G Z n , and where the initial data v[' F (0) = (^ r a (0)) QgZ n\.r } for 
1 < i < n and a € Z n \ {0} are given by n infinite vectors 

l ia) a eZ n \{0} ' 



vr' F (0) = = (^)L Z "Ucn . (94) 



Note that 

«io + r = 0, (95) 

such that we may cancel the zero modes. The justification for this is obtained 
for each iteration step m below. We may rewrite the equation (|93|) in the 
form 

<^ = A r > NS {v)v r ' F , (96) 
where v r ' F = (v^' F ,--- ,Vn F ^j . Furthermore A r,NS (v) is a nil 1 x nZ n - 

A r ' NS W = Uf s (v)) (97) 
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where for 1 < i,j < n the entry A 1 -'^ (v) is a Z" x Z n -matrix. We define 



A*" 8 (v) = [ ^ A[f S (v) vf , • • • , E (v) vf | , (98) 




where for all 1 < i < n 



T 



teZ™/ Ki<n 



The entries A^ a ^ (v) of A 7 "' (v) are determined by the equation in ([93 
of course. On the diagonal, i.e., for i = j we have the entries for a, (3 ^ 



Aaj/J ( V ) = <% E"=l " I W 1 )-^ Ej=l 



27TZQ!i l/^nl- ' 



(100) 



where for a = ft the terms of the form u^, « are zero (such that we do not 
need to exclude these terms explicitly). Furthermore, off-diagonal we have 
for % 7^ j the entries 

I 1 " *i)A^8 ( V ) = 2TOl{a^0} ^2 ■ ( 101 ) 

The idea for a global scheme is to determine v r,F = lim^oo v r ' m '^ for a 
simple control function and a certain iteration 



dv r ' m - F 
dt 



A NS ( v r,m-l\ v r,m,F ( 1Q2 ) 



starting with v r '° := h or with the infomation of a global solution of the mul- 
tivariate Burgers equation. The coefficients in (|102p will be time-dependent 
in general, which means that we shall have to develop a kind of Dyson for- 
malism for the equation. We shall define Trotter product formuls for certain 
linear infinite equations with time independent coefficients, and then we shall 
define a time discretization in order to apply these Trotter product formula. 
The time-dependent linear approximations are then solved by limits. 

In order to define a global solution scheme for this system of coupled 
infinite nonlinear ODEs we may start with the solution of the (viscous) 
multivariate Burgers equation (where we know that a unique global smooth 
solution exists), i.e., the solution of 
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where a G Z n , and where the initial data vf F (0) = (^(0)) z „ for 1 < 
i < n and a G Z n are given by n infinite vectors 



vf F (0) = hf = (^) T . (104) 

We know that the solution (v?) . is in /i s (Z n ) for all s G R for all 

time t > 0. However, without assuming knowledge about the multivariate 
Burgers equation we may also start with 



fit — " 2^ij=l \ J U ia 2-jj = \ Z^ 7 gZ" i ,L ]{a-l) U i~i 
+2 7 na;l{ c ^O} — 4to 2 ~ 

where a G Z n , and where the initial data v 0,i;, (0) = ^vJ'^(O), • • • , v^(0) 

are given by v°' F (0) = (f l ° a (0)) QgZ „ = (h ia ) a&n for 1 < i < n. In the 
equation (|105p and for v > the zero modes (\a\ =0) are the only modes 
where the damping (viscous) term 



(105) 



E 



Via (106) 



cancels. The same holds for the equation in (|82[) of course. Therefore we 
introduce the first approximation of a control function r^ : [0, oo) — > R 
for each 1 < i < n such the scheme for vp '^ = v^' F + eor° , 1 < i < n 
becomes a system of nonzero modes. At each stage we have to prove that 
this is possible, of course, i.e., that the solution of the linear approximative 

,0,F , „ J) _ („,0 \ _ f n r,0\ , „ 

with 



problem exists. Here, v-' + e o rf = (y i a ) aeZ n = (yl a J ^ + r? l0 is a vector 



I + r *o else. 

The control function is constructed in such a way that it cancels the 
zero modes of the original system at each stage of the construction. Hence 
it is constructed as a series (r™) m where the covergence asmjoo follows 
from properties of the controlled approximative solution functions ■v r ' m ' F . 
Note that it has to be shown that the control function is a) finite at each 
stage, and b) that it is finite in the limit of all stages. Note that the Leray 
projection term does not contribute to the zero modes in (|105p . and this 
may already indicates this finiteness. The next step is to analyze the scheme 
formally defined in the introduction starting from (|58|) to (|70|) . The effect 
of this formal scheme is that it is autonomous with respect to the nonzero 
modes, i.e., only the modes with \a\ ^ are dynamically active. That the 
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suppression of the zero modes via a control function is possible, i.e., that the 
controlled system is well-defined, is shown within the proof. Then we get 
into the heart of the proof. The iteration describes a set of analytic vectors 
(v^ F (t)) ew = U' m ' F (t)) with 

v r > m > F (t):=Tex V (A r m t)h(, (108) 

where T is a time order operator (as in Dyson's formalism), and the re- 
strictive dual function spaces together with the dissipative features of the 
operator A r m will make sure that at each time step m the approximation 
(|108|) (corresponding to a linear equation) really makes sense. Concerning 
dissipation features note that we may choose v > arbitrarily as is well- 
known and as we pointed out in the introduction. Note that we have a 
time-order operator T for all stages m > 1 since the the infinite matrices 
A r m depend on time t for m > 1. In order to make sense of matrix multi- 
plication of infinite unbounded matrices we consider the matrices involved 
in the computation of the first approximation v['°'' F (t) = Pq\^' F (t) (where 
Po is the related projection operator which eliminates the zero mode) as an 
instructive example. We may rewrite (|105|) in the form 

dt — u 2-^j=\\ I Via 2-rj=l 2^ 7 £Z" I ,L ]{a-l) U i~i 

+2ma-li , E ^=i E ^ ^(° t -7t)ft t ( a - 7) o (109) 



,V« 9<7r,7v-1r , 1 ^ gZ " 47r ^( a fc-7fc)fefc(^-7) „.0 



Hence this infinite linear system can be written in terms of a matrix Aq = Aq 
with n x n infinite matrix entries. More precisely, we have 

Ao=U j ) , (110) 

where for 1 < i, j < n we have 

A j = S ij D° + S ij C + L° ij (111) 
along with the Kronecker 5- function <5y, and with the infinite matrices 

„2 



D° := -^£?=l — 



£y - ^(27r«)ail{a^o} EE^f 



a„8€Z" 



22 



corresponding to the the Laplacian, the convection term, and the Leray 
projection terms respectively. Stricly speaking, we have defined a n x n 
matrix of infinite matrices where the Burgers equation terms are encoded 
on the n diagonal places and the linearized Leray projection terms exist also 
off diagonal. It is clear how the matrix multiplication may be defined in 
order to reformulate ()109p and we do not dwell on these trivial formalities. 
We note that 

ST ~ A o v > 

(113) 

v°> F (0) = h F 

along with the vectors = (vf, • • • , v^) T and = (h^, • • • , h F ^ T . Note 
that we have off-diagonal terms only because we consider global equations, 
i.e., equations, which correspond to linear partial integro-differential equa- 
tions. Note that the matrices D°,C ,L®j are unbounded even for regular 
data hi with fast decreasing modes (for the the matrix C° consider con- 
stant a — f3 and let /3j go to infinity for some j). The difficulty to handle 
unbounded infinite matrices in dual space becomes apparent. However, it 
is at this point that we can take advantage of the dissipative nature of the 
operator which is indeed the difference to the Euler equation. This is a 
major motivation for the dissipative Trotter product formula. Now, for two 
matrices M = (m a p) a a eZn and N = {n a p) a o eZn we may formally define 
the product P = {p ai ) an&n = MN via 

Pay = ^2 m <*P n Pi ( 114 ) 



for all a, 7 € Z n . There are natural spaces for which this definition makes 
sense (and which we introduce below). In order to apply this natural theory 
of infinite matrices developped below the next step is to observe that for 
k > matrices such as 



exp(^) (C°) fe = (exp (-vY.U ^) &=l 



'- h j(a-/3) 



\2mon l {c ^o} ElU4^f 



/ a,/3GZ" 

(115) 

have indeed bounded entries if the modes hip decreases sufficiently, and we 
shall see that with regular data hi related infinite matrix multiplications 
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are well-defined. Note that it is nevertheless diffcult to define the formal 
solution for v 0,i? , i.e., the expression 

exp ((SijD + SijC + L%)t) h F (116) 

directly. For this reason we shall consider projections of infinite vectors to 
finite vectors with modes of order less than I > 0, i.e., projections P v i which 
are define on infinite vectors (g Q )^ gZ n by 

P v l (9a)aeZ n := (9a)f a \<i (117) 

and projections P M i of infinite matrices ( m ai3) a ^ e x n 1° ^ n ^ e matrices with 

P M i (m a /5)ai3eZ" := ( m «/3)n<« (I 18 ) 

and prove Trotter product formulas for finite dissipative systems. This will 
lead us to the crucial observation then of a product formula of Trotter-type 
of the form 

lim ztoo exp ((<%P M! £>° + <%P M iC° + P M iL%)i} P v ih F 

= Km noo lim fctoo (exp {S^P^D^) exp ((c^P^C + P M iI%) ^ P Vl h F , 

(119) 

where (SijD°+5ijC +L^), i.e., the argument of the projection operator Pm ;) 
denote 'quadratic' matrices with (n x IT 1 ) rows and (n x Z n ) columns and 
h is some regular function. This means that for finite I > we can prove 
a Trotter type relation and in the limit the left side equals the solution 
v°' F of the equation above and is defined by the right side of (|119p . In 
the last step the regularity of the data comes into play. In this form 
this formula is useful only for the stage where the coefficient functions 
do not depend on time. For stages m > of the construction we have to 
take time dependence into account. However, we may define a Euler-type 
scheme and define substages which produce the approximations of stage 
m in the limit. The formula (|119p depends on an observation which uses 
the diagonal structure of D° (and therefore exp(Z) )). Here the matrix 

5ijC G +L^-—[ 5ijC° + L?- J , i.e., the deficiencies of symmetry in the matrices 

5ijC°+lJ-j, factorize with the correction term of an infinite analog of a special 
form of a CBH-type formula. It is an interesting fact that this deficiency is 
in a natural matrix space, we are going to define next. 

Concerning the natural matrix spaces, in order to make sense of formulas 
as in ()108p for s > n > 2 for a matrix M = (m Q . i a) Q ^ gZ n we say that 

M G (120) 
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if for all a, G Z n 

s i + | a %» < 121 > 

for some C > 0. For s > n and a vector it; = (w a ) ae z n G ^ s (Z) we define 
the multiplication of the infinite matrix M with w by Mw = (Mw a ) aeZn 
along with 

Mw a := m a pwp. (122) 
Indeed we shall observe that for r, s > n > 2 we have 

V < (123) 

^ \a - 8\ s 8 r ~ 1 + \a\ r+s - n y ' 

such that for s > n and M G we have indeed Mw G /i r (Z n ) if w G 
/i r (Z n ). This implies that for s > n, M G h s (Z n x Z n ) and io G h r {IT 1 ) 

M 1 w := Mw, M k+1 w = M w (124) 

is a well defined recursion for k > 1. Hence for a matrix M which is not 
time-dependent the analytic vector 

exp(M^ : =]T^^ (125) 

fc>0 

is well defined (even globally). For a time dependent matrix A = A(t) we 
formally define the time-ordered exponential 'a la Dyson' to be 

Texp(^t) := ^m=o mJ J[o t] dti--- dt m TA(t{) ■ ■ ■ A{t m )dti ■ ■ ■ dt m 

(126) 

:= £™=o lo dt i Jo 1 *2 • • • f- 1 ^i) • • • A(t m ). 

Note that T is the usual Dyson time-order operator which may be defined 
recursively using the Heavyside function. However, note that redefining 
with the Heavyside functions may make matters quite delicate if regularity 
with respect to time derivatives is considered. At each stage m > of the 
construction we shall define a scheme based on the Trotter product formula 
such that in the limit the linear approximation expressed formally by 

v r,m,F^ ._ Texp (j4 r^ fa F ^ 

gets a strict sense. This holds also for the uncontrolled linear approximation 

v m ' F (t) :=Texp(^ m t)h F , (128) 

where we shall see that these expressions can be well-defined for data h G 
h s (Z) for s > n and for t > 0. Note that proving the existence of a limit of 
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the the iteration with respect to m in a regular space also requires v > 0. 
In order to prove the existence of a limit there are basically two possibilities 
then. One is to prove the existence of a uniform bound 

< C (129) 

for some C > independent of m, and then proceed with a compactness 
arguments a la Rellich. A weaker form of (|129p without the time derivative 
and a strong spatial norm (large s) is another variation for this alternative 
since product formulas for Sobolev norms with s > \n and a priori esti- 
mates of Schauder type lead to an independent proof of the existence of a 
regular time derivative in the limit m f oo. Maybe the latter variation is 
the most simple one. An alternative is a contraction argument on a ball of 
an appropriate function space. Clearly, the radius of the ball will depend on 
the initial data, dimension and viscosity. This dependence can be encoded 
in a time weight of a time- weighted norm- as it is known from ODE theory 
for finite equations. Contraction arguments have the advantage that they 
lead to uniqueness results naturally. For (t — > w(t)) G C ([0, oo) x h s (Z n )), 
we define for some C > (depending only on v > 0, the dimension n > 0, 
and the initial data components hi G C°° (T n )) the norm 

\ w \hsc := SU P ex -P{-Ct)\w{t)\h°, (130) 

te[o,oo) 

and the norm 

I^I^C 1 : = SU P eM-Ct)(\w(t)\ h s + \D t w(t)\ h s), (131) 
te[o,oo) 

where D t w(t) := (jf^WcOagz™ denotes the vector of componentwise deriva- 
tives with respect to time t. Then a contraction property 

/ c- m.,F\ ( r,m,F /,\ r,m—l,F\ rion\ 

( SVi = w - v i ( 132 ) 

\ / l<i<n \ / l<i<n 

for all 1 < % < n with respect to both norms and for s>n + 2adl<i<n 
can be proved. Summarizing we have the following steps: 

i) in the frist step we do some matrix analysis. First, the multiplication 
of infinite matrices Aq and A m with infinite vectors h F is well defined. 
Second matrix multiplication in the matrix space is well-defined 
for s > n > 2 as is exp(M) for M G M*. Then we prove a dissipative 
Trotter product formula for finite systems and apply the result to the 
first infinite approximation system with solution v 0,i? at stage m = 
which is related to a linear integro-partial differential equation. 



sup 



■,m,F 



(t)\ 



h s + sup 



o 



r r,m,F 
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ii) In the second step we set up an Euler-type scheme based on the 
Trotter-product formula which shows that for some v > and s > 
n > 2 the exponential 

^ F (t):=(TeK V (A r m t)h F ). (133) 

is well-defined for all 1 < i < n, and in hf (Z n ) for all m, i.e., that the 
linearized equations for v[' m '^ which are equivalent to a linear partial 
integro-differential equation have a global solution. Here (.)j indicates 
that we project to the ith component of the n infinite entries of the 
solution vector at the first stage of construction. 

iii) for some v > as in step ii) the limit 

^ F (t):=(Te W (A^t)h F ). (134) 

exists, where 

Vi F (t) := lim mfoo v™' F (t) 

= (Texp (A^t) h F ). := lim mtoo (Texp (A r m t) h F ) i G hf (Z n ) . 

(135) 

This limit can be obtained by compactness arguments and by a con- 
traction results in time- weighted regular function spaces. The latter 
result leads to uniqueness of solutions in the time-weighted function 
space. 

Concerning the third step (iii) we note that for all t > and x S Tf we have 
v^it^x) := ^ v j£ hF (t) exp (27riax) (136) 

a&L n 

for all 1 < j < n. Hence, in classical Sobolev function spaces we have a 
compact sequence ( f^' m '' F (t, .) ) in higher order Sobolev spaces H r with 
r > s by the Rellich embedding theorem on compact manifolds for fixed 
t > 0. In this context recall that 

Theorem 2.1. For any g > s, q, s £ M and any compact manifold M the 
embedding 

j : if 9 (M) -> ^ s (M) (137) 

is compact. 

For any t > we have a limit .) £ H s (T n ) C C" 1 for s > m + f and 
the fact that the control function r is well-defined and continuous implies 

^.)=»J(t,.)- r (t)er (138) 
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for fixed t > and all m € N. Finally we verify that Vj is indeed a classi- 
cal solution of the original Navier-Stokes equation using the properties we 
proved for the vector v r,F in the dual formulation. One possibility to do this 
is to plug in the approximations of the solutions and estimate the remainder 
pointwise observing that it goes to zero in the strong norm |if|^ s • 



3 Proof of theorem 11.2 



We consider data hi E C°° (Tp) for 1 < i < n. The special domain of a re- 
torus has the advantage that we may represent approximations of a solution 
evaluated at a given time t > in the form of Fourier mode coefficients with 
respect to the basis 

i.e., with respect to an orthonormal basis of L 2 (T"). It is natural to start 
with an expansion of the solution v = (v±, ■ ■ ■ , v n ) T in the form 

Vi (t, x) = £ Q6Z n v ia exp (^p) , 

(140) 

P(t^)=Ea 6 ZnPae XP (2sp), 

where the modes Vi a , a € 7L n and p a , a € Z n depend on time. Here 
we have an index 1 < i < re which gives rise to some ambiguity with the 
complex number 2ni, but the latter i always occurs in the context of tt such 
that no confusion should arise. Plugging this ansatz into the Navier-Stokes 
equation system ([75]) we formally get re coupled infinite differential equations 
of ordinary type for the modes, i.e., for all 1 < i < n we have for all a € Z n 

dv ia _ v^n f _47ra| \ _ „ n „ 27ri7j 
dt — V L^i=\\ ~T~ j V ia 2^j=l Z^7GZ n I v j{a-~/) v ij 

V J (141) 

-2TTiaip a , 

where for each 1 < i < n and each «GZ™ 

Vi a ■ [0, oo) R (142) 
are some time-dependent a-modes of velocity, and 

p a : [0, oo) -> R (143) 

are some time-dependent a-modes of pressure to be determined in terms 
of the velocity modes Vi a . Here, (a — 7) denotes the subtraction between 
multiindices, i.e., 

(a - 7) = (ai - 71, • • • ,a n - 7 n ), (144) 
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where brackets are added for notational reasons in order to mark seperate 
multiindices. Next we may eliminate the pressure modes p a using Leray 
projection, which shows that the pressure p satisfies the Poisson equation 

= - ^2 v J,kVk,j, (145) 
i,k 

where we use the Einstein abbreviation for differentiation, i.e., 

dv ' 

Vj , k = ^ etc.. (146) 

The Poisson equation in (|145p is reexpressed by an infinite equation system 
for the a-modes of the form 

p J2 ~ i7r2(y i = jr A7llj ^ ak ~ lk)VjiVk[a - i) (147) 

i=l ^ j,k=l"/£Z n ^ 

with the formal solution (w.l.o.g. a/0 -see below) 

P*--t{ a ^0} ^ n _ x 47T 2 a 2 ' ^ } 

for a^O, which is indeed independent of the size of the torus I. This means 
that this term becomes large compared to the second order terms and the 
convection term for large tori (I > large while viscosity v > stays fixed). 
This may happen for approximations of Cauchy problems by equations for 
large tori for purpose of simulation. Note that 

' 1 if a ^ 0, 

!{«#>} := I (149) 
else. 

v 

Note that we put po — for a = in (|148|) . We are free to do so since 
(v,p + C) is a solution of (|76|) if (v,p) is. Plugging into (|14ip we get for all 

1 < i < n, and all aeZ™ 



dt ~ l^j=l v \ 12 J v ia 2_/j=l 2^ 7 eZ" l v j{a-i) v il 

(150) 

+ 27^1(^0} ~ Y\ n \ 4n 2 a 2 ' 

For 1 < i < n, this is a system of nonlinear ordinary differential equations 
of 'infinite dimension', i.e., for infinite modes Vi a and p a along with a S Z n . 
We emphasize that it may be rewritten in the basis 

cos (ax) , sin (ax) , a E Z n , (151) 
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with real modes at every step. Take (|150p as a shorthand notation for an 
equivalent equation with respect to a real basis as in (|15ip . This implies 
that we are constructing real solutions, although our notation is complex. 
For computations this real basis may be further reduced (cf. section on 
computational issues below). This means that we construct real solutions 
Vi a with Vi a (t) G R. Note that it suffices to prove that there is a regular 
solution to (|150|) . because this translates into a classical regular solution 
of the incompressible Navier-Stokes equation on the n-torus, and this im- 
plies the existence of a classical solution of the incompressible Navier-Stokes 
equation in its original formulation, i.e., without elimination of the pressure. 
We provided an argument of this well-known implication at the end of sec- 
tion 1. This is not the controlled scheme which we considered in the first 
section. However, this scheme is identical with the controlled scheme at the 
first stage m = 0, so we start with some general considerations which apply 
to this first stage first and introduce the control function later. So let us 
look at the simple scheme in more detail now. The first approximation is 
the system 

dt ~ 2-^=1 I "—p- i u ia z^j=iZ^ 7 eZ' 1 i '^(0-7)^7 

(152) 



Let 



Ej,fc=i £ 7 ez" 47r 7j(Qfc-7ft>° 7 frfc( a - 7 ) 
fi,F _ ( 0,F 0,F\ 0,F ._ / \T 



K'\---^ F ), vr ( 153 ) 

F 

Formally, we consider v- ' as an infinite vector, where the upper script T in 
the componentwise description indicates transposition. Then equation (|152p 
may be formally rewritten in the form 

^r+£4-vr> ( i54 ) 



dt 



with the infinite matrix A l n = ( a*° fl ) and the entries 

a a/3 - 6 *pV I - 2^=1 -p - I ~ 2-,j=l Z 1" (155 J 

and where 

L ijaP = l{am 27Tla i ^2 ( 156 ) 

for all 1 < i,j < n denote coupling terms related to the Leray projection 
term. As we observed in the preceding section you may write (|154p in the 
form 

<iv ' F . „ P , , 

— = • < 157 > 
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where 



4) V 



0,F ._ 



((^■a^)+L°^)^ , (158) 



\ Ve{l,-,n},/3GZ" / aeZ" y 



Ki<n 



and j4o = (^o ) can ^ e considered as a quadratic matrix with n 2 infinite 
matrix entries 

J$ = Al for i = j, (159) 

and 

Al j = L% for i ? j. (160) 

Note that these matrices which determine the first linear approximation 
matrix Aq are not time-dependent. The modes a^g,^, or at least one set 
of these modes, have to decay appropriately as |a|,|/3| t 00 i n order that 
the definition in (|158p makes sense, i.e., leads to finite results in appropriate 
norms which show that the infinite set of modes in Aqv 0,f belong to a 
regular function in classical space. Since the matrix Aq has constant entries, 
the formal solution of ()154p is 

v°' F = exp{A t)h F . (161) 

In order to make sense of this formula we shall use the dissipative feature 
on the diagonal terms and a Trotter-type product formula (otherwise we 
are in trouble because the modulus of entries increases with the order of 
the modes). We solve the equation (|264p by approximations via systems of 
finite modes, i.e., via 

P vl v°' F; * = exp {PmM) hF (162) 

where P v i and P M i denote projections of vectors and matrices to finite cut- 
off vectors and finite cut-off matrices of modes of order less or equal to I. 
Note that we used the notation 

Pviv /'* (163) 

with a star upperscript since the projection of the solution of the infinite 
system is not equal to the solution of the projected system in general. The 
latter is an approximation of the former which becomes equal in the limit. 
Indeed we use 

v°' F = linn too P v iv° ,F '* = linn too lim fctoo ^exp (P M i (<%£> ) |) x 

(164) 

\ k 

x exp 
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Note that exp ((<%£>°) £) exp ^(SijC + L? J J |J e M*. Here the sym- 
bol (.),y indicates projection onto the infinite 7* n x Z"-block at the ith row 
and the jth column of the matrix Aq. This means that the right side is 
well defined due to the matrix exp ((5ijD°) £) which has the effect of an 
multiplication of each row by a function which is exponentially decreasing 
with respect to time and with respect to the order of the modes. The more 
delicate thing is to prove that the Trotter-type approximation converges to 
the (approximative) solution at each stage m of the construction. We come 
back to this later in the proof of the dissipative Trotter product formula. 
Let us start with the description of the other stages m > of the construc- 
tion first. At stage m > 1 having computed v™ -1 '^ = (v^ -1 ) QgZ n the rath 
approximation is computed via the system 



dt ~ l^j=l \ ~! r ~ J ia ^j=l ^ 



76Z™ I v j( a -y) v iy 



(165) 



with the same initial data, of course. We have 

vl^ == («K&3P. , (166) 
and the equation (j!65|) may be formally rewritten in the form 

with the infinite matrix A m = ( a 1-7 ^ 1 ) and the (time-dependent!) en- 

tries A m = = (a*f) , and along with A m = D° + Cg + L£ for all 

1 < i < n, and Am = L™ for i ^ j, 1 < i,j < n, where for all i 7^ j we have 

rm _o_ ; - n . ^=1 4 ^( a *-^)"M°-^) (1681 



2 "7j „m-l ,,m 



The matrix C™ related to the convetion term — Y^=i S 7 ez™ 
is defined analogously as in the first stage but with coefficients instead 
of hi a . Again, the equation (|166p makes sense for regular data and in- 
ductively assumed regular coefficient functions because polynomially 
decay of the modes hi a and Vi a compensate that encoded quadratic growth 
due to derivatives up to second order. For a solution however we have to 
make sense of exponential functions with exponent defined in terms 
of modes v^~ l and applied to data h F . Again we shall use a dissipative 
Trotter product formula in order to have appropriate decay as \a\, \/3\ f 00 
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in order for the solution v^. However, this time we have time-dependent 
coefficients and we need a subscheme at each stage m of the construction in 
order to deal with time dependent coefficients. The formal solution of (|167p 
is 

v m ' F = Texp(A m t)h F , (169) 

where Texp(.) is a Dyson-type time-order operator T defined above. Note 
that for all m > 1 the functions v™ 1 ^, 1 < i < n represent formal solutions of 
partial integro-differential equation if we rewrite them in the original space. 
Again in order to make sense of them we shall use the dissipative feature and 
a Trotter-type product formula at each substage which is a natural time- 
discretization. It seems that even at this stage we need viscosity v > in 
order to obtain solutions for these linear equation in this dual context. This 
assumption is also needed when we consider the limit m | oo. We shall 
estimate at each stage m of the construction 

v m ' F = Texp(A m t)h F (170) 

based on the Trotter product formula where we set up an Euler-type scheme 
in order to deal with time-dependent infinite matrices in the limit of sub- 
stages at each stage m. 

Next we go into the details of this plan. Let us consider some linear alge- 
bra of time-independent infinite matrices with fast decaying entries (which 
can be applied directly at the stage m = 0). The matrices of the scheme 
considered above are (n x Z n ) x (n x Z n )-matrices, but - for the sake of 
simplicity of notation - we shall consider some matrix-algebra for Z n x Z n - 
matrices. The considerations can be easily adapted to the formally more 
complicated case. First let D = (d a /3) aj3eZ „ and E = (d a p) a p eZ „ be two 
infinite matrices, and define (formally) 

D ■ E = (/ aj9 ) aj8ez „ , (171) 

where 

fa/3 = ^2 (172) 
■y£Z n 

Next we define a space of matrices such that (|17ip makes sense. For s € R 
we define 

< := = {d a p) af3&n | d a p e C & 3C> : \d a p\ < 1 + ■ 

(173) 

In the following we consider rather regular spaces where s > 2 + n. Some 
results can be optimized with respect to regularity, but our purpose here is 
full regularity in the end. Next we have 

Lemma 3.1. Let D £ and E € MJ for some s, r > n + 2. Then 

D ■ E £ M r n +S ~ n (174) 
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Proof. For some c > we have for a//J 



I fa/3 1 — I X^eZ" d a ^e y i3\ 

c 

^7^{o,/3} |a-7| s |7-/3| r 



^ c + X^{a,/3} | a — y|« | 7 - C /3|'- ( 175 ) 



< C+ cC 



|a-/9| r+s - n 

The latter inequality is easily obtained b comparison with the integral 

dy 



»\B al) \a-y\ r \y-P\ f 



(176) 



where B a p is the union of balls of radius 1/2 around a and around /3. Partial 
integration of these intergals in polar coordinate form in different cases leads 
to the conclusion. We observe here that there is some advantage here in the 
analysis compared to the analysis in classical space because we can avoid 
the analysis of singularities in discrete space. □ 

For r = s > n this behavior allows us to define iterations of matrix 
multiplications recursively according to the matrix rules defined in (|17ip 
and (|172p . i.e., we may define by recursion 

A = /, where I = (5 a p 
A 1 = A, (177) 
A k+1 = A ■ A k . 

In the matrix space we may define exponentials. For our problem this 
space is too 'narrow' to apply this space directly. However, it is useful to 
note that we have exponentials in this space. 

Corollary 3.2. Let D £ for some s > n. Then 

eX p(£>) = ^— eM* (178) 

fc=0 

Well what we need is an estimate for the approximative solutions. In 
this context we note 

Corollary 3.3. Let D G M* for some s > 2n and let h F G h s (Z n ). Then 

exp(L>)h F G h s (Z n ) . (179) 
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Proof. Let 

F = exp(D), where F = (/ Q/3 ) Q(3eZ n, (180) 

and let 

9a = ^2 f a P h P- ( 181 ) 



Then for some C, c > 

In? — 

7GZ" /3eZ"\{a,0} 

□ 



a.1 = I £ < £ i^-^i^ ^tth 7 - (182) 



We cannot apply the preceding lemmas directly due to the fact that the 
diagonal matrix with entries — udij Y!k=i 17~ a k * s no * bounded. Neither is 
the the matrix related to the convection term. However the multiplication 
of the disspative exponential with iterative multiplications of the convection 
term matrix stay in a regular matrix space such that a multiplication with 
regular data of polynomial decay lead to regular results. The dissipative 
term - the smoothing effect of which is obvious in classical space - makes the 
difference. At first glance, iterations of the matrix lead to matrices which 
live in weaker and weaker spaces, and we really need the dissipative fea- 
ture, i.e., the minus signs in the diagonal in mathematical terms, in order 
to detect the smoothing effect in the exponential form. The irregularity re- 
lated to a positive sign reminds us of the fact that heat equations cannot be 
solved backwards in general. However, due to its diagonal structure and its 
negative sign we can prove BCH-type formula for infinite matrices. The fol- 
lowing has some similarity with Kato's results for semigroups of dissipative 
operators (cf. [H]). However Kato's results (cf. [H]) usually require that 
the domain of the dissipative operator includes the domain of the second 
operator summand, and this is not true in our formulation for the incom- 
pressible Navier-Stokes equation, because the diagonal operator related to 
the Laplacian, i.e., the diagonal terms (— 5 a pv Y^,i=l a i ) a / 9ez™\{o}' i ncrease 
quadratically with the order of the modes a. So it seems that the cannot be 
applied directly. Anyway we have 

Lemma 3.4. Let gf be finite vectors of modes of order less or equal to I > 0. 
Then for some finite vector ff = (f a )i a \<.i with finite entries f l a , and in the 
situation of Lemma \3. 7| we have 



^exp (CH) gf - exp [{A 1 + B l )t) gf ) J < \f l a t% (183) 

where (.) denotes the projection to the ath component of an infinite vector. 

This is what we need in the following but let us have a closer look at the 
Trotter product formula. First we define 
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Definition 3.5. A diagonal matrix ($ a pd a p) a g eZ „ is called strictly dissipa- 
tive of order m > if d aa < — c\a\ m for all a € Z n and for some constant 
c > 0. It is called dissipative if it is dissipative of some order m. 

In order to state the Trotter-product type result we introduce some no- 
tation. 

Definition 3.6. For all m > 1 let 7 m = (7™, ■ ■ ■ 7^) denote multiindices 
with m components and with nonnegative integer entries 7™ for m > 1 and 
1 < i < m. For each m > 1 we denote the set of m-tuples with nonnegative 
entries by Nq\ Let B = (b a p) a ^ e im denote a quadratic matrix, and let 
B T = (bpa) a g^n be its transposed, and E be some other infinite matrix of 
the same type. For m > 1 and m-tuples j m = (7™, ■ ■ ■ )7™) with nonneg- 
ative entries jj 1 , 1 < j < n, we introduce some abbreviations for certain 
iterations of Lie brackets operations of matrices. These are iterations of the 
matrix AB := B — B T and either the matrix B or the matrix B T in arbi- 
trary order. This gives different expressions dependent on the matrix with 
which we start, and we define J 7 m (starting with [AB, B]) and 1^ (starting 
with [AB,B T ] I? accordingly. First we define Jy (AB, B) = J ( i)(A5, B) 
(starting with [A£?,i? T ] for 71 > recursively). Let 

[E, B] T = EB — B T E, (184) 

which is a Lie-bracket type operation with the transposed. For 7^ = define 

I bl) [AB,B]:=AB, (185) 

and for 7]- > define 

J (7 i) [AB, B] := [J 71 _! [AB, B] , B\ T . (186) 
Having defined E m. — l and if 7™ = then 

7 7 m [AS, fl] = I 7 m-i [A5, B] + 7 7 m-i [A5, B T ] (187) 
Finally, if 7™ > 0, then define 

E r [AB, B] = [l jm -i [AB, B],B] T . (188) 
Similarly, for 7J = define 

and for 7J > define 



lT[AB,B):=AB, (189) 



[AB, B] := [I 71 _! [A5, £ T ] ,B) T . (190) 
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Having denned I~ m — 1 and if 7™ = then 

I$n [AB, B] = lT r -, [AB, B] + 7^_! [AS, 5 T ] (191) 
Finally, if 7™ > 0, then define 

I^m [AB, B] = [l^-i [AS, 5] , S T ] . (192) 

Next we prove a special CBH- formula for finite matrices. We do not need 
the full force of the lemma [3771 below for our purpose, but it has some interest 
of its own. The reader who is interested only in the global existence proof 
may skip it and consider the simplified alternative considerations in order 
to see that how a Trotter product result can be applied. The results may 
be generalized but our main purpose in this article is to define a converging 
algorithm for the incompressible Navier-Stokes equation which provides also 
a constructive approach to global existence. We show 

Lemma 3.7. Define the set of finite modes 

ZP := {a G Z n ||a| < I}. (193) 

Let A 1 be the cut-off of order I of the dissipative diagonal matrix of order 2 
related to the Laplacian and let B l = (b a /3) a p € %n be the cut-off of some other 

matrix. Next for an arbitrary finite quadratic matrix N = (n a p), a i let 

exp m (TV) = exp (AO - HT — . (194) 

Vfc=o ' / 

Then the relation 

exp(A l ) exp(S') = exp(C z ), (195) 
holds where C l is of the form 



(196) 



C l = A 1 + B l + \ exp(2A0A^ + Yni>\ e xp m (A 1 ) x 
x (cpmlpm [AB l ,B l ] T + C T m lJ m [AB',J3 l ] T ) 
for some constants cpm , cj m such that the series 

£ cpmt m (197) 



m>l 



and the series 



E c V tm ( 198 ) 



m>l 



converges absolutely for all t > 0. 
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Proof. For simplicity of notation we suppress the upper script I in the fol- 
lowing. All matrices A, B,C, • • • are finite multiindexed matrices of modes 
of order leass or equal to I. The matrix C = Yli^i Q 1S formally determined 
via power series 



C(x) = ^2dx\ C'{x) = J^iCjx* -1 
i=l i=l 

by the relation 

OO 

]T R k [C'(x), C(x)} =A + B + Y,R h [AB]^ 



(199) 



(200) 



k=0 k>l 

where for matrices E, F [E, F] = EF — FE denotes the Lie bracket and 



R l [E, F] = [E,F], R k+1 [E,F] = \R k [E,F],F 



(201) 



recursively (Lie bracket iteration on the right). Comparing the terms of 
different orders one gets successively for the first C-terms 

Ci = A + B, C 2 = \ [A, B] , 
C S = ±[A, [A,B}} + ±[[A,B},B], 
C 4 = ±[A,[IA,B],B}} + ±[[A,IA,B}},B] 
C 5 = [A [[A, [A, B}} ,B}} + ^ [A, [[[A, B] , B}} , B] 
- 3±o [A, [[[A, B],B], B)\ - gig [[A, [A, [A, B])} , B) 



(202) 



UA, [A, [A, [A, B]]]] 



720 



720 



[[[A,B},B],B],B],--- 



Iterated Lie brackets simplify if A is a diagonal matrix. First we define left 
Lie-bracket iteratons, i.e., 



L l [E,F] = [E,F], L k+l [E,F] 



E 



L k [E,F] 



(203) 



recursively. Next we study the effect of alternative applications of iterations 
of the left and right Lie-bracket operation for this specific A. We have 



[A, B] = AAB, 



(204) 



with AB = B-B T , and 



L fc [A,B] = A k 2 k ~ 1 AB, 



(205) 
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Next we have 

R k [A, B] = AR^ 1 [AB, B] T (206) 

Next 

LR k [A, B] = A 2 R k ~ 1 ([AB, B] T + [AB, B T ] T ) (207) 

Induction leads to the observation that other summands than A + B of the 
series for C can be written in terms of expressions of the form 

[A) k Ip m [AB,B] T (208) 

and in terms of expressions of the form 

(A) k lJ m [AB,B] T (209) 

For each order p = k + m we have a factor with < TP summands. This 
leads to global convergence of the coefficients cpm , c^ m . □ 

We continue to describe consequences in a framework which is very close 
to the requirements of the systems related to the incompressible Navier- 
Stokes equation. As a simple consequence of the preceding lemma we have 

Lemma 3.8. Let gf be a finite vector of modes of order less or equal to 
I > . In the situation of Lemma \3.1\ we have 



lim (exp (a 1 ~\ exp (^)) 9i = e*P {{A 1 + B l )t) gf . (210) 

It is rather straightforward to reformulate the latter result to equations 
with matrices of type P m iAq, P M iA r Q , i.e. finite approximations of order I 
of the nil 1 x nZ"-matrices Aq, Aq. We have 

Corollary 3.9. Let B^ 1 = ( P M i (c? + L?) ) and 



ij 



B r b ' l '° = (P* (Cf + such that A l = D° b > l +B /, A r ' 1 = D°/+B r b ' 1 . 

Then the Trotter product formula 

i^M^lM^))'^-*)- (211) 

and the Trotter product formula 

BS h ( D >''1) exp " exp K<) (212) 

hold. 
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Next note that iterations of B?' 1 of order k > 1 have at most linear 
growth for constellations of multiindexes where a — 7 is constant (with the 
order of the modes as I f 00). Let is have a closer look ar this. We have 



(213) 



We consider the entries for the square. We have 



Z^i=i —r~ h 3(a-P) + 27rzail{a^o} 4^2^ J (214) 

2^j=i — r~ ft i(/3-7) + ^A 1 ^} EILi 471-2/^ 

Expanding the latter product (|214p and inspectin the four terms of the 
expansion we observe that two of these four terms are entries of matrices 
in the regular matrix space. Only the mixed products of convection terms 
entries and Leray projection terms entries lead to expressions which are a 
little less regular. We define an appropriate matrix space 

M*' lin : = 

(215) 

{D = (d a p) a p ezn d l af3 G C & 3C> V/ : \d l a p\ < C +ffgy } • 
By induction we have 

Lemma 3.10. Given 1 < i, j < n for all k > 1 we have 

B b i, a ,) k ^M^ n (216) 

Hence the statement of lemma 13.101 is true also for the exponential of 
£{V.. We note 

Lemma 3.11. If D is a strictly dissipative diagonal matrix of order 2, then 
for geven 1 < i,j <n we have 

exp(D)exp(s 6 %.^) eM*. (217) 

In order to establish Trotter product representations for our iterative 
linear approximations v r,m,F of a controlled Navier-Stokes equation we con- 
sider a matrix space for sequences of finite matrices {D l ) l where each D l 
has modes of order less or equal to I. We define 

M*' fin := 

(218) 



C 



^ G C&3O0VZ:|^|< ITra 
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Here for each I > 1 the set denotes the set of modes of order less or equal 
to I. 

It follows that for all t, k and 

e X p(^0exp(^^))%M^ (219) 

for s > n, i.e., the finite approximations of the Trotter product formula are 
regular matrices indeed. Hence 



feJtl exp v° ; i)^^ '*)) £M » (220) 

Now consider finite approximations of the problem (|157p to the cut-offs 
of this problem to a set of problems of finite modes with modes of order less 
or equal to I, i.e., 



For each I the solution 

v°< F = exp hf (222) 



is globally well-defined via the dissipative Trotter product formula. Our 
infinite linear algebra lemmas above imply that the sequence (v^i'j is a 
Cauchy sequence for s > n if h F E h s (Z n ) for all set Hence we have 

Lemma 3.12. We consider the torus of length I = 1 w.l.o.g.. Let h F E 

h s (Z n ) for allseR. 

v^=(ex P (A t)h^). 

/ k x (223) 

= hmjtoo (hm fetoo (exp (l^f) exp (B b ' l j;)) hf J E h" (Z") 



whenever hf E h s (Z n ) /or s > n for 1 < i < n. 

For the same reason 

Corollary 3.13. Consider the same situation as in the preceding lemma. 

v['°< F = (exp(^)h^ 

/ k \ ( 224 ) 

= lim, too ^lim fctoo (exp(^ W |)exp(^ W |)) hf J E h s (Z n ) 

whenever hf E /if (Z n ) /or s > n for I < i < n. 
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Proof. Let Aq = P m iAq denote the projection of the matrix Aq to the finite 
matrix of modes less or equal to modes of order I. For finite vectors v^f ,F 
with 

r,0,F r,0,F 
v i = j mv M ( 225 ) 

we have 




(226) 



The right side of (|226|) is a Cauchy sequence with index I in M*. □ 

The stage m = is special as the matrix ij£' '' does not depend on 
time. For this reason we get similar Trotter formulas for time derivatives 
using the damping of the dissipative factor exp (D r b ' 0tl \ We shall use first 
order time derivatives formulas for linear problems with time-independent 
coefficients later when we approximate time dependent linear problems via 
time discretizations. We have 

Corollary 3.14. Recall that Aq 1 = P m iAq denotes the projection of the 
matrix A r to the finite matrix of modes less or equal to modes of order I . 

K'°' F = l(expW)^) j 

= lirm too (V hm fctoo (exp (l^f) exp (fif ^))* hf \ G h s (Z n ) 

(227) 

whenever hf € hf (Z n ) for s > n for 1 < i < n. 

Next at stage m > 1 we cannot apply the results above 
to define 

(v m > F ) =Te W (A m t)h F , 

or in order to define 

( v r> m > F )=Texp(A r m t)h F . (229) 

The main difference of the stages m > 1 to the stage m = is the time 
dependence of the coefficients formally expressed by the operators T (|228[) 
and (|229p above. The second difference is that we need to control the zero 
modes if we want to establish a global scheme. For the latter reason, next 
we consider the controlled scheme (the considerations apply to the uncon- 
trolled scheme as well for one iteration step). In each iteration step we 
construct the solution of an infinite linear ODE in dual space which corres- 
peonds to a linear partial integro-differential equation in the original space. 



directly in order 
(228) 
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For each iteration step we need some subiterations in order to deal with the 
time-dependence of the coefficients. In order to apply our observations con- 
cerning linear algebra of infinite systems and the Trotter product formula, 
we consider time-discretizations. The time dependent formulas in (|228p and 
(|229p have rigorous definition via double limits (with respect to time and 
with respect to modes) of Trotter product formulas for finite systems. There 
are basically two possibilities to define a time-discretized scheme based on 
this form of a Trotter product formula. One possibility is to consider the 
successive linerarized global problems at each stage m of the construction, 
where the matrices A r m are known in terms of the entries of the infinite vector 
v r,m-i,F w liicla contains information known from the previous step. Accord- 
ing to the dissipative Trotter product formula we expect a time-discretization 
error of order 0(h) where h is an upper bound of the time step sizes. The 
other possibility is to apply the dissipative Trotter product formula locally 
and establish a time local limit lim^oo \ r ' m ^ F at each time step. We shall 
consider the first possibility. 

In any case, we define a time-discretized subscheme at each stage m > 1 
and use the Trotter product formula for dissipative operators above locally in 
time in order to show that we get a converging subscheme. The convergence 
can be based on compactness arguments, and it can be based on global con- 
traction with respect to a time-weighted norm introduced above. Another 
possibility is to establish a time-local contraction at each stage m > 1, i.e. 
in order to establish a solution v r > m ' F atstage m. The time- local contrac- 
tion can be iterated due to the semi-group property of the operator such 
that the subscheme becomes global in time. At each iteration step m the 
time-dsicretized scheme for each global linear equation works for the uncon- 
trolled scheme and the controlled scheme if we know that the data v m_1 ' 
(resp. v r,m-1,ir ) computed at the preceding step are bounded and regular, 
i.e., v r ' m ~ 1 ' F (t) £ h s (Z n ) for s > n + 2 with n > 3. Both subschemes have 
a limit with and without control function. It is in the limit with respect 
to the stage m that the control function r becomes useful. It also becomes 
useful for designing algorithms. Next we write down the subscheme for the 
uncontrolled subscheme and the controlled subschemes. Later we observe 
that the control function is globally bounded in order to prove that there 
is a limit as m f oo. For each step m of the construction, however, the 
arguments for the controlled and the uncontrolled scheme are on the same 
footing. We describe the global linearized scheme, i.e., the scheme based 
on global linear equations which are equivalent to partial integro-differential 
equations in the original space. The procedure is defined recursively. For 
t > arbitrary we define a scheme which converges on [0, t] to the approx- 
imative solution v r ' m '^. Having defined v m ~ 1,F at stage m > 1 we define 
a series v r ' m ' F,p , p > 1, where p is a natural number index which refers to 
a global time discretization tq, p > 1 and where q £ {1, • • • 2 P } along with 
tq — = 2~ p for all p > 1, and to := 0, such that we get a contractive 
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scheme with respect to a certain time- weighted Sobolev norm. Next for 
given p > 1 and g £ {!,-•■ , 2 P } we define the equation for ■v r ' m > F >P>Q on the 



interval 



L q-l> L Q 



and where the initial data are given by 



r r,m,F,p,q-l [ g> ^ 



Here for q = 1 we have the initial data 

v r,m,F,p,0 (^ 



h 



(230) 



(231) 



Note that the vector h r ' equals the inital data h , but without the zero 
modes, i.e., 



where for 1 < i < n we have 



T 



ti> F 



(h ia ) 



a£Z n \{0} 



(232) 



(233) 



Next we define for each p > 1 the sequence of local equations for v r,m,F ' p,q . 
This leads to a global sequence (v r,m,ir ' p ) >l defined on [0, T] for arbitrary 
T > 0. The next goal is to obtain a contraction result with respect to 
the iteration number p in order to establish the existence of global regular 
solutions v r ' m ' F for the approximative problems 



dv r,m,F 



dt 



j^r,NS /y.r,m—l\ ^ r r,m,F 



L) v rv 



(234) 



where v r,m,F 



r,m.F 



r.m.F 



T 



and where for each m > 1 the initial 



data are given by h r,F . Here we look at A r ' N defined above as an operator 
such that A r ' NS (v r,m_1 ) is obtained by applying this operator to the argu- 



ment v r,m instead of v r . We shall give the details for the iterations steps 
p,q > 1. First we have to obtain \r r > m > F > p > q for each p > 1 and q G {1, • • • 2 P }. 

we have the local Cauchy 

j^r,NS ( v r,m-l^ v r,m,F,p,q^ (235) 



On the time interval 
problem 



f P f P 



2p ^ 5 2P 



dt 



where v r ' m ' F < p < q = ^ v r 1 ' m ' F ' p ' q ; . . . ; v ^ m ' F ' p ' q J and where for each m > 1 the 

initial data are given by v r ' m ' F ' p ' q ~ l f • Next we explicitly describe the 
matrix A r ' NS (y r > m - 1 >P> q ) which is a n (Z n \ {0}) x nZ n \ {0}-matrix with 



j^r,NS / r,m— 1 



) = « 5 (v— 



l<i,j'<n 



(236) 
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where for 1 < i,j < n the entry Alf S (v^" 1 " 1 ) is a Z n \{0} x Z™\{0}-matrix. 
We have 

j^r,NS ^ v r,m— 1^ ^r,m,F,p,q _ 

A r ' NS ^ v r,m-l^ v r,m,F,p,q _ _ _ j^r,NS (^ r r.m-l\ ^r,m,F,p,q\ T 



jj=l^ij V v J v l ' > Z-,j=l ^ry 

where for all 1 < i < n 

\ \ r 

r,WS / rm _i\ r,m,F,p,q\ 



En v-> at, ISO I 

7 = 1 Z^fleZ™ ^iai'fl l V 



ioj/3 V v ; V jf3 



(237) 



(238) 



The entries -A^L-g (v r,m x ) are determined by the equation as follows. On 
the diagonal, i.e., for i = j and for a, (3 ^ we have the entries 

°V A iaj/3 V I - °n 2^j=l v I J °*j2^j=l I V j(a-/3) 

(239) 

2vrzai l{ a ^o} 47rQ ^ 1 , 

where for a = j3 the terms of the form are zero (such that we do not 

need to exclude these terms explicitly). Furthermore, off-diagonal we have 
for % 7^ j the entries 

Ylk^^Pjiak - /3 fc )r r '' 



(1 - SirfA^p (v) = 2vr?aa {Q ^ 0} v „ , — 2 ■ ( 240 ) 



The idea for a global scheme is to determine v r,F = lim^oo v r ' m '^ for a 
simple control function and a certain iteration 

dv r, m ,F _ NS , r m -l\ „r,m,F 



dt 



A NS f v r,m-l\ v r,m,F^ ( 2 41) 



starting with v r, ° := h. We shall use the abbreviation 

A r m := A NS (^ m ~ 1 ) . (242) 

The proof of global regular existence can be obtained if a uniform up- 
per bound for the sequence iy r,m,F ) m> i can De found. An alternative is a 
contraction result. Contraction results have the advantage that they lead 
to uniqueness results with respect to the related Banach space. Next we 
define a Banach space in order to get a contraction result for the subscheme 
(v r ' m ' F ' p ) >r Note that for each p > 1 the problem for v r ' m ' F > p on [0,T] 



is 



defined by 2 P recursively defined subproblems for v r ' m,F,p, ' } for 1 < q < 2 P 
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which are defined on the interval t p where the data for the problem 

for q > 2 are defined by the final data of the subproblem for v r,m,F ' p ~ q ~ 1 
evaluated at t p v 

Next consider the function space 

B n,s ;= | t ^ u F {t) = ( U F . . . >U £)T ^ > ^ ^ £ ^ {%n ^ {Q}) j 

(243) 

Note that we excluded the zero modes bcause this Banach space is designed 
for the controlled equations. For u F E B n,s define 



u 



T,exp 

h a ,C 



sup Vexp(-Ct) uf(t) 
te[o,T) i=1 



(244) 



In the following we abbriviate 



A r 



A 



NS (, r r,m-l,F 



(245) 



Especially the evaluation of the right side of (|245|) at time t is denoted by 
A r m (t). Recall that 

z n,0 ._ z n \ | Q | ^ 



For s > let us assume that 

t v p ' m - 1 ' Ji '(t) G C fc ([0,T],/i s (Z n '°) .) 



(247) 



Here C fc ([0,T],/i s (Z n '°)) is the function space of n x (Z n, °) vectors which 
have /c-times differentiable component functions t — > w^ m_1 (i) for 1 < i < n 
and a G Z n \ {0}. Note that there is no need here to be very restrictive 
with respect to the degree of the Sobolev norm s. Each component function 
satisfies a Taylor formula. Hence we have the Taylor formula 



_r,m— 1,F 



{t + h) = Y.«< P <k-iD P t v 



r,m—l,F 



(t)hP 



+T^Jo 1 U-*) 
for te[0,T-h] and h > 0. Here 



k— 1 j~jkyr,m— 1,F 



(248) 



r.m— 1,F 



W : = ( ^ 



r,m— 1 



In the following we assume that that h^' , v, 

n + 2 and 1 < i < n, because this property is inherited for v^'" l,r as we 
go from stage m — 1 to stage m. Note that the matrix- valued function 



l<i<n,a£Z n '° 
r,m—l,F 



(249) 



€ /i s (Z n -°) for 



r,m,F 
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t -> A r m (t + h) - A r m {t) = A NS (v r ' m ~ 1 ' F (t + h)) - A NS (v r ' m - 1 - p (t)) is 
Lipschitz, i.e., for some finite Lipschitz constant L > we have 



(250) 



< L|v r ' m - 1 '- F (t + ft) - v r ' m - 1 ' F (t)| fc .. 
Using the infinite linear algebra lemmas above we observe 



Lemma 3.15. Let s > n + 2 and T > be given. For k = 2 assume that 
v r,m-i,F £ S re g U i ar as in [2^1). Then for some finite C > 



su P«e[o,T] exp(-Cu) 



r,m,F,p-l 



a) 



< 



h» 



2''- 



(251) 



Proof. For notational reasons the size of the torus is assumed to be one. We 
remark that 



and accordingly 



^m(*2(g-l)) ~~ ^m(^g-l) 



(252) 



(253) 



For some vector C f : € /i s (Z n ) for s > n + 2 we postulate the difference 



r r,m,F,p,q-l ( ,p 



2(9-1) 



r,m,F,p-l,q-l ( .p-1 
'9-1 



■,h,p 
J q-1, 



and we consider some properties which the vector Cq' p inherits from C™^ 
Next at stage p > 1 consider the initial data v r ' m,F (t^Z\) of the problem 



(254) 



ih,p 



at substep 1 < q < 2 P 1 . We have t € 
makes sense to consider the the subintervals 



l 2(q-l)' l 2q 



b q-H L Q 



f P f P 



and 



where it 



f p f p 

L 2q-li L 2q 



We may consider t 6 



f p f p 

L 2q-li L 2q 



w.l.o.g. because the following estimate 



simplifies t G 



t v f p 

l 2q-2i l 2q-\ 



For t £ 



f P f P 
L 2q-1' L 2q 



we have 



r,m,F,p,2q 



(*) " vl 



r,m,F,p—l,q 



< 



r,m,F,p,2q 



r,m,F,p,2q i,p 



r,m,F,p—l,q 
i 



r,m,F,p—l,q i,p 



(tlq-S) 



+ 



r,m,F,p,2q-l,.p 
[Z 2q 



(Vi)-v^ 



r,m,F,p—l,q/,p 



(Vi) 



(255) 
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Since _ 1 \ = t P q -\ anc ^ with (|254[) above for the last term in (|256p we have 



+ 



V; 



r,m,F,p,2q—l / ,p 



(&-l) ~ < 



r,m,F,p—l,q /,p 



(Vi) 



(exp (4(^1)) " Vi))) vr ' m ' F ' P ' 29 (^( 9 -i)) 

(exp (^(t£) - i^)) v^- 1 -^-^(^: i 1 )) i 



(256) 



< 



exp ^m(*2(g-l) 

exp {%) C% 



for some finite C > 0, which depends only on data known at stage m — 1. 
Furthermore, for the first term on the right side of (|256p we may use the 
rough estimate 



r,m,F,p,2q 



h s 



< 



exp (^(^_!) (t - vr^^J - v^^^,) 



exp 



f^m(^2g-l) 



f - £ 



2g-l 



r,m,F,p—l,q /,p 



( Vi) - < 



r,m,F,p—l,q /,p 



The right side of (|257|) we observe with (|256p 



(257) 



exp ( ( t - t* ) ) v^' 2 ^*) - vj 



r,m,F,p—l,q /,p 



(Vi) 



< 2 



r,m,F,p,2q r.p 



(Vi)-V 

exp (f ) 



r,m,F,p—X,q/.p 



(Vl) 



(258) 

As for each p > 1 the entries of the sequence Cq' p are in O (/i 2 ) where 
/i denotes the maximal time step size (which is 2~ p with our choice) the 
diiference to be estimated is in 0(h) and we are done. 

□ 

The preceding lemma shows that we have a Cauchy sequence 

( vr ' m ' F ' P ) P >i ( 2 59) 
with respect to a regular (time weighted) norm and with a limit v r,m ' F with 



v[' m ' F (i) G h s (Z n >°) 



(260) 
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for all t G [0,T]. Since we have a dissipative term (damping exponential) in 
the Trotter product formula, similar observations can be made for the time 
derivative sequence 

d r ,m,F, P \ , . 

Note, however, that s > n + 2 > 5 above may be arbitrary, and this can 
be exploited in order to prove regularity with respect to time t for each 
^r,m,F y j a defying equation of the latter function. We even do not 
need estimates for products of functions in Sobolev spaces which may be 
borrowed from classical Sobolev space analysis. Using the lemma above for 
large s > we may instead use the regularity implied by infinite matrix 
products as pointed our above. As a consequence of the preceding lemma 
we note 

Lemma 3.16. For all m > 1 and s > n + 2 > 5 the function 

(v wF ) , = (Texp (A r m t) h r ' F ) i G h\ (Z*°) , (262) 

is well-defined, whenever hf G hf (Z n ). 

The same holds for the uncontrolled approximations, of course. We note 
Corollary 3.17. For all m> 1 and s > n + 2 the function 

(^ F ) i = (Texp (Ant) h F ) t G hf (Z n ) (263) 

is well-defined, whenever hf G (Z n ). 

However, it is essential to get a uniformly bounded sequence 

K' m ' F L e N = ( Tex P hf e ^ (^)) me N • ( 264 ) 

for some ^ > 0. We remarked in the introduction that we can even choose 
v > 0. Indeed, we observed that we can choose it arbitrarily large (as is 
also well-known). However this was not needed so far and we shall not need 
it later on. It is just an useful observation in order check algorithms in the 
most simple situation via equivalent formulations with rigorous damping. 
At this point it is useful to consider the controlled sequence ( v r -' m ' F ) 

V /m€N, l<i<n 

Recall that the functions v r ' m '^ are designed such that the zero modes are 
zero i^q" 1 = 0. The control function is just defined this way. Next we show 

that a uniformly bounded controlled sequence ( v^ ,m,F ) implies 

V /men, l<i<n 

unformly boundedness of the uncontrolled sequence (v"^) . In 

V /meH, l<i<n 

order to observe this we go back to (|2ip . At stage m > 1 it is assumed the 
v^q 1 1 = 0. The controlled approximating eqation at stage m is obtained 
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from (|2ip by elimination of the zero modes. We have for 1 < i < n and 
a/0 the equation 

ifcJJ™ v- ,n | 47ra| \ r m v - ^ 27ri7j r,m— 1 r,m 

— 3T~ ~~ 2^j=l 17 I i 5- J V ia ~ 2^j=l 2^ 7 gZ' 1 \{0,q} I V j(a--y) V i-y 

M=l E 76 Z"\{o, a } 47T7j(Qfc-7fc)^ 7 « fc(ct _ 7) 
+^*«i 1 {o#0} E Li ' 

(265) 

We considered the controlled equation systems as autonomous systems of 
non-zero modes. However, in order to compare the controlled system with 
the original one we may define 

dv iQ m — V n V 271-ryj r,m-X r,m f 2 fifil 

Note that vj' m— ■v r i ' m,F e h s (Z n '°) for s > n > 3 implies that the right 
side of ()266|) is bounded by constant finite C > 0. Hence we have 

Lemma 3.18. If the sequence (v r ' m ' F ^ has an upper bound c > with 
respect to the \.\ s = Yli=i \-\h"-norm, we have for some finite C > 

\r (t)\ < exp(Ct) (267) 

for all t > 0. 

It remains to show that for some finite C > such that for all 1 < i < n 
and m > we have 

K' m > F (t,.)\ hf <C (268) 

Based on the arguments so far there are several ways to get an uniform 
bound for the sequence (v r,m ) One is via contraction results on cer- 
tain balls in appropriate function spaces. The radius of such a ball clearly 
depends on the size of the initial data and on the size of the horizon. How- 
ever it is sufficient that for each dual Sobolev norm index s > and for 
each data size \h F \ and each horizon size T > we find a contraction re- 
suit on an appropriae ball for a related time weighted function space. Let's 
look at the details. For arbitrary T > consider two smooth vector- valued 
functions on the n-torus, i.e., functions of the form 

f,g£ [C°°([0,T],T n )] n . (269) 

Consider the equations 

dv^ =A r,NS^ v r,f,F f ( 270 ) 
along with V'f' F (0) = h r,F , and 

*!%£ = A r > NS ( g )v r >°> F , (271) 
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along with v r ' 9 ' F (0) = h r ' F . Here we denote v r 'f' F = (v^' F , ■ ■ ■ ,Vn^' F ^j 

and similarly for the function v T,s,J? . As in or notation above the matrix 
A r,NS ( f ) is a nZ n x n z n -matrix with 

\r,NS s e \ _ ( A r,NS 



A r ^(f) = (^f a (f) (272) 
where for 1 < i, j < n the entry A^ 3 (f) is a Z n x Z n -matrix. We define 

In n \ 

A r,NS (f) v r,f,F = £ ^5 (f } yf (f ) , (273) 

\i=i i=i / 

where for all 1 < i < n 



aeZ™/ Ki<r 



The entries (f ) of (v) are determined as follows. On the diago- 

nal, i.e., for i = j we have the entries for a, (3 ^ 

<Miaj/j ( f ) = <% Di=i v [ — v^j-^j Di=i —r L fj(c t -fi) 

(275) 

-l-Ojj Z7Tiaj l{ a ^0} f« 47rQ ,y j 

Z-^i— 1 j 

where for a = f3 the terms of the form fu a -p\ are zero (such that we do not 
need to exclude these terms explicitly). Furthermore, off-diagonal we have 
for i 7^ j the entries 

(i - (f) = 2.»,i 1 ^ ol ^i±|g0M^. (276) 

The definition of A r ' NS (g) is analogous. Next for functions u r ' F = \vQ F , • • • , 
and for s > n + 2 consider the norm 



u 



r,F\ Y^n \,J' > F 



Er=il u i L.- (277) 



Consider a ball of radius 2|h r ' | around the origin, i.e., consider the ball 

I s 

For s > n + 2 and for data h r ' € /i s (Z™ \ {0}) the considerations above 
shows that the Cauchy problem 

d^' F =A r > NS (f)v r >f> F , (279) 
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along with v r '^ ,i? (0) = h r,F has a regular solution v r, f' F with v r "' (i) 6 
/i s (Z™ \ {0}). Moreprecisely, for s = n + 2>5we have 



W rJ ' F \T„ < 2\h r ' F \ T. 



(280) 



Next for linear operator 



(g))u 



iaj/8 



.ryZVS 



iajfS 



r,F 



(281) 



are Lipschitz with some Lipschitz constant L on this ball. Note that we have 



A r,NS (f) _ A r,NS , x 



(282) 



9 • -, Efc = i47r/3 j (a fc -/3 fc )(/ fc(Q _^ ) -g fc(c »_ / 3 ) ) 



Furthermore, for g r,F 6 ^ 2 |h rF | 
constant L is chosen such that 



we may assume w.l.o.g. that the Lipschitz 



sup g ,,F eB | (SijA^f (g)) u^. < L 



u 



r,F I 



(283) 



where the weaker norm on the right side of (|283p is due to the fact that we 
have to compensate the first term on the right side of 



47TOj 



(284) 



, _- i H2=i in Pj( a k-l3k)gk( a -i3) 
+2Tria>il {am £ ^ 4 ^ a 2 • 



Lemma 3.19. Let T > be arbitrary, and let P' , g^* g £? 

s > n + 3. For C > 31 we have 



2\h r < F \ 



for 



y r,f,F _ v r,g,F 



T, exp i 
< - 

h°,C 2 



g 



T, exp 
h",C 



(285) 



For C > 6L we have 



v r,f,F _ V ,g,F 



T,exp,i 1 
< - 
h",C 2 



g 



T, exp, i 
h",C 



(286) 
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Proof. For each t £ [0, T] we have 

V r </> F (t) - V r '9' F {t) 



< 



< 



+ 



Jt A r,NS ( f ) (u)V r J> F (U)du - jj A r ' NS (g) (u)v r '9' F (u)du 

(j£ A r ' NS (f ) (u) - J * A r ' NS (g) («)) v^(u)d« 

jj A r ' Ar ' s (g) («) (v^ F (u)ds - v r '»' F (u)) rf n 

< LTsup ue[0ir] f F (u) - g F (u) 

+LT sup ug [ 0T ] v r, f' F (u)ds — v r ' 9 ' F (u 
It follows that 

V r ^ F (t) - V r '9' F (t) 



< L 



r,exp 



f - g / Q exp(Cu)d' 



h a ,C 
v rJ,F _ v r,g,F 



T,exp . 

L exp(Cu)du 
h s ,c JU 



< 



Lcxp(Ct) 
C 



f F "g F 



T,exp 

h s ,C 



t C su Pu>0 



v r >f' F (u)ds - v r '9'- F (.) 



Since 



v r,f,F _ v^.f.-F 1 



T,exp 



< 



it follows that 



v rJ,F _ v r,g,F 



r,exp 



T,exp 

h a ~ 2 ,C 

T,exp 



h s ,C 



< 



1 L 



f F (u)-g F (u) 



T,exp 

h s ,C 



(287) 



(289) 



(290) 



For C = 3L the result follows. The reasoning for the stronger norm is 
similar. □ 
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It is clear that this contraction result leads to global existence: for each 
s > n + 2 and T > we find a constant C > such that 



r,m,F 



T,exp 

< C (291) 

h°,C 



Uniqueness follows by standard arguments. Uniform upper bounds for the 
approximative (controlled) solutions v r ' m,F lead to existence via compact- 
ness as well. Note that the infinite vectors v^ ,m ' F (t) = (t)) aeZ n,o are in 
1-1 correspondence with classical functions 

v? m (t,x)= <r^p(^), (292) 

where up m £ iP (T n ) for s > n + 2. Recall that 

Theorem 3.20. For r > s and for any compact Riemann manifold M (and 
especially for M = Tf) we have a compact embedding 

e : H r (M) -> H s (M) (293) 

This means that (^[' m ) mg pj has a convergent subsequence in H r (Tf) for 
r > s which corresponds to converging subsequence in the corresponding 
Sobolev space of infinite vectors of modes. Hence, passing to an appropriate 
subsequence v[' m ,F (t) of vl' m ' F (t) we have a limit 

vf(t) = lim v^Ct) € (Z») (294) 

Tn'yoo 

for r < s (Rellich embedding). Since s is arbitrary this limit exists in h r (Z n ) 

for all r£R. Hence, for all 1 < i < n we have a family ( v, r,m ' F ) which 

V /m'eN 

satisfies 

rff^ m / 47ra? \ r m ' _ n _ 27U7j r,m' — 1 r,m' 



En a i \ r.m —1 r.rn 

+2TOl{ c ^0} E? W • 

(295) 

If we can prove that the limit of infinite vectors ( — ^ — (t) ] is contin- 
uous in the sense that 



dt 



V / a eZ« (296) 

:= { (9a) | E ae z« /« «cp G C (Z»)} , 
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then we obtain a classical solution. Inspecting the terms on the right side 
the assumption that s > 2 + n and n > 2 is more than sufficient in order to 
get 

Aira 2 

—p 1 - ) <<f ' I e /i s " 2 (Z n ) C /i n (Z n ) , (297) 



such that with this assumption we may ensure that I — 1| — I (i) con- 

V / mGN 

verges in C (Z n ) C /i r (Z n ) with r > if we can control the expressions for 
the convection term and for the Leray projection term appropriately. How- 
ever, this is easily done with the help of the infinite linear algebra results 
above. We stick to s > n + 2 and n > 2. First for the convection term for 
each 1 < i < n and all a £ Z n we consider 

-E E ( 29s ) 

j=i 7 ez"\{«} 

We observe |7jf[^, m |/ l «-i(z») ^ C f° r some constant C > independent of 
m, hence 

E ') u» g /i2s " 1 " n * c ( 2 ") 

for some C > independent of m such that the limit is in h 2 (Z n ). Hence 
(|298p and the limit for m' f oo is in h? (Z n ). Similarly, the Leray projection 
term 

En ,i / \ r.m' — 1 r.m' 

3 -,fc=i E 7£Z " 47r 7j K ~ TfcHr ^ (a - 7 ) 

is bounded by a product of two infinite vectors which have some uniform 
bound C > in h s_1 (Z n ), such that the Leray projection term is safely 
in /i 2 ( s ~ 1 )~ n (Z n ) C h 2 (Z n ) where we did not even take the |a| 2 in the 
denominator corresponding to the Laplacian kernel into account. We have 
shown 

Lemma 3.21. For s > n + 2 ami n > 2, and for the same v > as above 
there is a C > such that for all 1 < i < n and m' > we have 



h 



dt 



uniformly for t > 0, and 



±vf{t) = lim Ur'' F {t) € V (Z») C C (Z») C ^ (Z») . (302) 
dt m'too at 
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We conclude 



Theorem 3.22. The function 

v[' F (i) = lim vp ro '' F (t), 1 < i < n (303) 

m'foo 

satisfies the infinite nonlinear ODE equivalent to the controlled incompress- 
ible Navier-Stokes equation on the n-torus in a classical sense. Moreover, 
since the argument above can be repeated with arbitrary large s > we have 
that for all 1 < i < n the infinite vector v^ F (t) and its time derivative are 
in h s (Z n ). Higher order time derivatives also exist in a classical sense by 
an analogous argument for derivatives of the Navier-Stokes equation. 

Finally we have 

Theorem 3.23. Let h { G C°° (¥"). For each v > and I > and for all 

1 < i < n and all t > 

Vi (t,.)eC^(Tf). (304) 

and 

v[(t)eh s (T?). (305) 

for arbitrary s£R. 

Proof. We have vf (t) = v\' F \t) - r(t) £ h s (Z n ) for all s G M, because 
v i' F (t) G h s (Z n ) for all s € R, and r(t) is a constant. The second part 
follows from Corollary above. Given s > we can differentiate 

i(t,.) = ^ 7 exp f^P^J), (306) 



up to order m, where m is the largest integer less than s, and get a Fourier 

/ 



series which converges in L? (T™). Hence 



Vi(t, .) G H s (Tf) , (307) 

for all t > 0. Since, this is true for all for all s > we the first statement of 
this lemma is true. □ 

Now for a fixed / > and any hf G h s (Z n ) for s > + 2 we have 

^v7^(t)e^(Z») (308) 



for all m > from the unifrom bound 



v. 



(*)Uf < C. (309) 
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4 A converging algorithm 

The infinite ODE-system can be approximated by finite ODE systems of 
with modes of order less or equal to some positive integer k in the sense of 
a projection on the modes of order \a\ < k. At a fixed time t > this is 
a projection from the intersection D s ^h s (Z n ) of dual Sobolev spaces into 
the space of (even globally) analytic functions (finite Fourier series). The 
Trotter product formula for dissipative operators considered above holds 
for the system of finite modes, of course. More importantly, the resulting 
schemes based on this Trotter product formula approximate the correspond- 
ing scheme of infinite modes if the limit A; y oo is considered. In this section 
we define this approximating scheme of finite modes. It is not difficult to 
show that the error of this scheme converges to zero as k f oo. Detailed 
error estimates which relate the maximal mode of the finite system and di- 
mension n and viscosity v to the error in h s norm are of interest in order to 
design algorithms. This deserves a closer investigation which will be done 
elsewhere. In this section we define algorithms of of different approxima- 
tion order in time via finite-mode approximations of the non-linear infinite 
ODEs which are equivalent to the incompressible Navier-Stokes equation. 
Furthermore some observations for the choice of parameters of the size of 
the domain / > and the viscosity v > are made, and we observe re- 
ductions to a cosinus basis (symmetric data). We consider the controlled 
scheme which solves an equivalent equation for the sake of simplicity. From 
our description in the introduction it is clear how to compute an approxi- 
mate solution of the incompressible Navier-Stokes equation from these data. 
Recall the representation of the controlled Navier-Stokes equation in terms 
of infinite matrices. It makes sense to formulate an 'algorithm' first for the 
infinite nonlinear ODE system. Accordingly, in the following the use of the 
word 'compute' related to substeps in an infinite scheme is not meant in a 
strict sense. It may be defined in some more strict sense b use of transfinite 
Turing machines, but what we have in mind is the finite approxmations of 
an infinite object such that schemes can be implemented eventually. This is 
not an algorithm in a strict sense (even the description is not finite), but the 
projection to a set of finite modes leads to the algorithm immediately once 
it is described in the infinite set-up. Consider the limit v r,F = (w[ Q ,) Q6Z „ of 
the construction of the last section, i.e., the global solution of the controlled 
infinite controlled ODE system, as given. Then we may write this function 
formally in terms of the Dyson formalism as 

v^ F {t) = T exp (AH) hf G hf (Z n ) , (310) 

where 

Texp(A r t) := E"=o mT I [0 1] dt i ■ ■ ■ dt k TA r (h) ■ ■ ■ A r (t k )dh ■ ■ ■ dt k 

(311) 

:= IXo Jo *i Jo 1 *a ■ ■ ■ Jo*" 1 *"(*i) • • • 
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where 

A-fc) := ,. r „ (312) 

and for all 1 < i, j < n we have 

A r ' ij = (a% j ) (313) 



where 



3' = ^ S^v - E - =1 ^ - ELi ) + ^ (314) 



^ = ^ 1 2 T ^ (315) 



along with 

Efe=l 4 7r ftj( a fc -Pk )Vk(a-P) 

Er=i47r 2 ^ 

Note that the modes t>j Q are time-dependent - for this reason the Dyson time- 
order operator appears in this formal representation. In order to construct a 
computable approximation of this formula we need a stable dissipation term 
for the second order part of the operator. We can achieve this if we may 
use the Trotter product formula above, because the factor exp {{5ijD r )) is 
indeed in the regular matrix space and it is a damping factor for the 
controlled system. In order to apply the Trotter product formula to the 
controlled Navier-Stokes system we may consider a time discretization and 
apply this formal time step by time step. Starting with v^' to := hi let us 
assume that we have defined a scheme for = To < t\ < ■ ■ ■ < i*. and 
computed the modes ' for < I < i& Then in order to define the scheme 
recursively on the interval [t^, ifc+i] we first consider the natural splitting of 
the operator and define 

A r (t J ) = (5 lJ D) + {B r ^(t k )), (316) 
where B r, ^(tk) = \b r ' t g' lJ ) along with 



^=(-EU 2m " jV r- p) ) +r%, (317) 



and 



Lg = 27TU* - EILi47f2a , (318) 

The most simple infinite scheme we could have in mind is then a scheme 
approximating v r ' F (t) = {vl a {t)) a ^ n for arbitrary given t > in m steps 
then is 

(exp (± (<%£>)) exp (f (B^(t m ))))x 

(319) 

.•.(ex P (^(%D))ex P (|(5^(i ))))h. 
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This is indeed a time-dependent Trotter-product type approximation of or- 
der o (-M J, as m t 00, where o denotes the small Landau o. Higher order 
approximations can be achieved taking account of some correction terms in 
the Trotter product formula. Indeed, consider 'truncations of order q } of the 
term defined in (|196p above. For each t r , < r < m consider 

C*ttr) = (SijD) + (BW(t r )) + 

Hp=l pi Sj>l, I3 1 £N[ , m'+\/3 l \=p, Km'+l X (320) 

xl^, [A (B r ^(t r )),(B r ^(t r ))] T . 

Acording to the Trotter product formula higher order achemes can be ob- 
tained if there is a correction term E such that the replacement of (B r ' lJ (t m )) 
by (B r M (t m )) + E cancels the correction of order q at each time step which 
may be defined by 

C 9 (tr)corr = C«(t r ) - ((%£>) + (B r ^(t r )) 

= Sp=l pi J2l>l, /3 l eN[ , m>+\p'\=p, Km'+l x (321) 

xlpj [A {B^{t r )),(B^{t r ))] T . 

The explicit analysis of these schemes has a right in its own an will be 
considered elsewhere. For simulations we have to make a cut-off leaving 
only finitely many modes, of course. The formulas established remain true 
if we consider natural projections to systems of finite modes. We define 



(322) 



(exp (±P Ml (5 i3 D)) exp (|P M! {B^(t m )))) x 

• • • (exp (±P Ml (S^D)) exp (£P M * {B™(t )))) P vl h F , 
where the operator P v i is defined by 

P v ih F = (P lv ihf , ■■■ , P nv ih F ) (323) 
such that for 1 < i < n and hi = (hi a ) ae %n we define 

P iv i(v a ) a&n = U) , (324) 

and for infinite matrices M = (M^ ) = ( ( m\-, ) ) we define 

P Ml M = (P ijM iMV) , (325) 
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where 



) 



\a\,\/3\<k 



(326) 



In both cases (vector- and matrix-projection) we understand that the order 
of the modes is preserved of course. Next it is a consequence of our existence 
constructive proof that 

Theorem 4.1. Let hi S C°°(Ti, and let T > be a time horizon for the 
incompressible Navier-Stokes equation problem on the n-torus. Let s > n + 2 
Then the finite mode scheme defined in \322\i converges for each < t < T to 
the solution of incompressible Navier-Stokes equation in its infinite nonlinear 
ODE representation for the Fourier modes as m, I f oo. 

Next recall that the parameter v > can be chosen arbitrarily. A large 
parameter v > increases damping and makes the computation more stable. 
However, in order to approximate a Cauchy problem (the initial value prob- 
lem on the whole space) we need large I (length of the torus) . Note that the 
damping terms in the computation scheme are of order fj, the convection 
terms are of order j and the Leray projection terms are independent of the 
length of the torus. They become dominant if the size I of the n-torus is 
large. Prom our constructive existence proof we have 

Corollary 4.2. The Cauchy problem for the incompressible Navier Stokes 
equation may be approximated by a computational approximation of order k 
on the n-torus Ti for I large enough where in a bi-parameter transformation 
of the original problem v > should be chosen 



such that the damping term is not dominated by the Leray projection terms. 
Furthermore due to the quadratic growth of the moduli of diagonal ( damping ) 
terms with the order of the modes compared to linear growth of the convection 
term and the Leray projection term with the order of the modes the scheme 
remains stable for larger maximal order k if it is stable for lower maximal 
order k. Furthermore the choice 



for s > 2 + n leads to solutions which are uniformly bounded with respect to 
time. 

For numerical and computational purposes it is useful to represent ini- 
tial data in cos and sin terms (otherwise computational errors may lead to 
nonzero imaginary parts of the computational approximations). We have 

Lemma 4.3. Functions on the n-torus Tf be represented by symmetric data, 
i.e. with the basis of even functions 



v>l 



2 



(327) 



v>2\h\ 2 s 



(328) 




(329) 
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The reason for the statement of lemma 14.31 simply is that you may con- 
sider a L 2 - function / on the cube [— l,l] n (arbitrary prescribed on [0, l] n ) 
with periodic boundary conditions with respect to a general basis for L 2 (Tf) 
such that 



and then you observe that /(•••, Xj. •••) = /(••• , — X{. • • • ) imply that the 
sin-terms implicit in the representation (|330p cancel. Hence on the n-torus 
which is built from the cube [0, l] n we may consider symmetric data without 
loss of generalization. 

5 The concept of turbulence 

In the context of the results above let us make some final remarks concerning 
the concept of turbulence. There are several authors (cf. [2] and [3] for 
example), who emphasized that the dynamical properties of the Navier- 
Stokes equation such as the existence of strange attr actors tec. rather than 
the question of global smooth existence may be important for the concept of 
turbulence. We share this view, and we want to emphasize from our point of 
view, why we consider the infinite dynamical systems of velocity modes to be 
the interesting object in order to start the study of turbulence. This concept 
is indeed difficult, and in a perspective of modelling it is always worth to 
follow this difficulty up to its origins on a logical or at least very elementary 
level. According to classical (pre-Fregean) logic, a concept is a list of notions, 
each of which is a list of notions and so on. This is not a definition but 
something to start with (the definition has some circularity as there is no 
sharp distinction between a concept and a notion, but this may be part of an 
inherent difficulty to define the concept of concept). According to Leibniz 
a concept is clear (german: 'klar') if enough notions are known in order 
decide the subsumption of a given object under a concept, and a concept 
is conspicuous (german: 'deutlich') if there is a complete list of notions of 
that concept. All this is relative to a 'cognoscens' of course. Since Leibniz 
expected that even for empirical concepts like 'gold' the list of notions may 
be infinite he expected a conspicuous cognition to be accessible only to an 
infinite mind and not to human beings. Anyway, according to this concept 
of concept we may say that our concept of the concept of turbulence may not 
even be clear. However, there is some agreement that some specific notions 
belong to the notions of 'turbulence'. Some aspects of fluid dynamics may 
be better discussed with respect to more advanced modelling. For example 
for a realistic discussion of a 'whirl' we may better define a Navier-Stokes 
equation with a free boundary in one (upper) direction and a gravitational 
force in the opposite direction (pointing downwards). This free boundary 
may be analyzed by front fixing on a fixed half space as we did it in the 




(330) 
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case of an American derivative free boundary problem. Additional effects 
are created by boundary conditions (the shape of a river bed, for example). 
Indeed without boundary conditions dissipative features will dominate in 
the end. First we observe that it is indeed essential to study the dynamics 
of the modes. The Navier-Stokes equation is a model of classical physics 
(|v| < c, where c > denotes the speed of light. Therefore physics should 
be invariant with respect to the Galilei transformation. If v = (vi)J <i<n 
is a solution of the Navier Stokes equation then for a particle which is at 
x o = O^oi) " " " j x 0n) at time to we have the trajectory x(t), t>to determined 
by the equation 

x (t) = v, i.e., Xi(t) = xoi + / Vi(s)ds, 1 < i < n, (331) 

Jto 

which means that for the modes we have 

Xai + / Vi a (s)ds, 1 < % < n, (332) 
Jt 

Hence, rotationality, periodic behavior, or irregular dynamic behaviour due 
to a superposition of different modes i>j Q which are 'out of phase' translates 
to an analogous dynamic behavior for the trajectories with a superposition 
of a uniform translative movement (which disappears in a moved laboratory 
and is physically irrelevant therefore, up to Galiliei transformation dynamic 
behavior of the modes translates into equivalent. The dynamics of the modes 
is an infinite ODE which may be studied via bifurcation theory. Since we 
proved the polynomial decay of modes a natural question is wether a center 
manifold theorem and the existence of Hopf bifurcations may be proved. 
The quadratic terms of the modes in the ODE of the incompressible Navier 
Stokes equation in the dual formulation seem to imply this. Irregular behav- 
ior may also be due to several Hopfbifurcations with different periodicity. 
Bifurcation analysis may lead then to a proof of structural instability, sensi- 
tive behaviour and chaos (cf. [1]). Global sensitive behavior and chaos may 
also be proved via generalisations or Sarkovskii's theorem or via intersection 
theory of algebraic varieties as proposed in [TJ] and p3] ■ Studying the effects 
of boundary conditions and free boundary conditions for turbulent behavior 
may also be crucial. For the Cauchy problem dissipative effects will take 
over in the long run. For a 'river' modelled by fixed boundary conditions 
for the second and third velocity components v 2 and V3 (river bank) and a 
free boundary condition (surface) and a fixed boundary condition (bottom) 
for the first velocity component v\, and a gravitational force in the direction 
related to the first velocity component, complex dynamical behavior can be 
due to the boundary effects. It is natural to study such a boundary model 
by the front fixing method considered in [5]. It is likely that the scheme 
considered in [7] and [H] can be applied in order to obtain global existence 
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results. 
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